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along  the  satellite  track  or  directly  beneath  the  satellite.  This  study 
considers  techniques  by  which  a  more  general  areal  survey  can  be  made  so  that 
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20.  ABSTRACT  (Continued) 
techniques  is  considered. 

Topside  sounders  currently  offer  the  best  areal  survey  capability  with 
the  ISIS  satellites.  A  spacecraft  version  of  the  "Digisonde"  sounder  offers 
enhanced  measurements. 

Laser  radars  or  Ildars  have  substantial  potential  for  ionospheric  mea¬ 
surements  but  have  not  yet  been  developed  for  spacecraft  use  or  for  operation 
in  the  vacuum  ultraviolet  region  of  the  spectrum  where  ionic  species  such  as 
0'*'  must  be  detected.  It  is  anticipated  that  such  developments  are  likely. 

Thomson  scatter  radars  are  useful  in  ground-based  installations  but 
would  be  cumbersome  on  a  satellite.  They  are  huge  and  require  high  power  to 
operate  effectively. 

Various  Faraday  and  phase  methods  have  been  used  over  the  years  but 
yield  only  the  column  electron  density  and  generally  require  a  co-operating 
ground  station.  They  are  useful  in  conjuction  with  other  measurements  but 
probably  not  alone. 

Numerous  optical  or  chemiluminescent  methods  are  possible.  They  often 
use  relatively  simple  equipment  and,  hence,  would  be  easy  to  include  with 
other  detectors  on  a  satellite  system,  v’ 

In  situ  methods  offer  no  possiblity  of  an  areal  survey  but  can  be  use¬ 
ful  where  detailed  composition  and  high  accuracy  are  important. 

Bremsstrahlung  is  an  indirect  method  potentially  useful  in  the  high 
latitude  regions  only.  Be  measuring  the  bremsstrahlung  x-ray  flux  associated 
with  a  diffuse  aurora  it  may  be  possible  to  determine  the  flux  of  precipitat¬ 
ing  electrons.  This  can,  in  turn,  be  related  to  the  ionospheric  electron 
density  in  regions  where  the  precipitating  flux  provides  the  dominant 
ionization. 

Because  no  single  method  offers  a  complete  picture  of  the  ionosphere, 
a  hybrid  method  involving  more  than  one  technique  is  clearly  preferable. 
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SECTION  1 

SUMMARY  AND  CONCLUSIONS 


A  variety  of  techniques  for  the  remote  sensing  of  ionospheric 
electron  density  from  a  satellite  platform  have  been  studied.  Because 
most  current  satellite  techniques  are  limited  to  determining  ionospheric 
characteristics  along  the  satellite  track  or  directly  beneath  the 
satellite,  this  study  considered  techniques  by  which  a  more  general  areal 
survey  can  be  made  so  that  the  electron  density  in  three  dimensions  can  be 
determined.  Table  1  provides  a  concise  summary  of  the  results. 

Topside  sounders  offer  the  best  areal  survey  capability  at  the 
present  time.  Both  of  the  ISIS  spacecraft,  launched  into  polar  orbits 
over  ten  years  ago,  examine  most  of  the  topside  of  the  earth's  ionosphere 
twice  each  day.  But  they  only  observe  during  a  small  interval  of  local 
time  and  without  the  sophisticated  sounding  techniques  now  possible.  Use 
of  a  spacecraft  version  of  the  "Digisonde"  sounder  on  several  polar  orbit¬ 
ing  satellites  would  considerably  enhance  topside  sounding.  This  instru¬ 
ment  is  capable  of  extracting  most  of  the  information  available  in  the 
radio  wave  reflected  from  the  ionosphere.  Phase  measurements  increase  the 
group  height  resolution  from  15  km  for  ISIS  to  0.1  km  for  the  Digisonde. 
Doppler  measurements  permit  determination  of  the  motion  of  individual 
ionospheric  layers. 

Laser  radars  (lidars)  have  considerable  potential  for  iono¬ 
spheric  studies  but  have  not  been  developed  to  the  point  where  they  could 
be  deployed  on  a  satellite  in  the  immediate  future.  These  instruments  are 
not  limited  to  viewing  the  subsatellite  track  beneath  the  spacecraft,  and, 
hence,  come  the  closest  to  providing  true  areal  coverage.  To  be  successful 
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Table  1.  Summary  of  techniques  for  measuring  ionospheric  electron 
density  from  a  spacecraft. 


Method 

Technique 

Electron 

Density 

Availability 

1—  0)^-0 
<o  >  fO 

Q)  Cl. 

J-  3  ro 
c  to  O 

Resolution 

Comment 
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Sensing 
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V 

C 

F 

G 

Best 

Current 

b)  Digisonde 

V 
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Technique 
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E 
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Future 

Technique 
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V 

FU 

VG 

F 
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4)  Faraday  and 

CL 

C 

F 

VG 

Needs 

Phase  Methods 

(V) 

Ground 
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Passive 

5)  Chemiluminescent 

CL 

C 

E 

VG 

Many 

Sensing 

Methods 

(V) 

Possi- 

bilities 

6)  In  Situ  Methods 

V 

C 

P 

E 

7)  Bremsstrahlung 

V 

NF 

G 

F 

Indirect 

Method 

8)  Hybrid  Methods 

V/CL 

FU 

E 

E 

Many 

Possi- 
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Codes:  V  =  Volume 

CL  =  Column 

(V)  =  Some  Volume  Information 
C  =  Current  Technology 
NF  =  Near  Future 
FU  *  Future 


E  =  Excellent 
VG  =  Very  Good 
G  =  Good 
F  =  Fair 
P  *  Poor 


in  measuring  F-region  electron  density,  lidars  would  have  to  detect  O'*', 
because  this  is  the  dominant  ionic  species  in  the  F-region.  It  would  be 
necessary  to  use  0^  resonant  or  ground  state  transitions  which  occur  in 
the  vacuum  ultraviolet.  Because  a  tunable  lidar  does  not  yet  exist  in  this 
region,  space-borne  lidars  for  the  study  of  the  F-region  are  hypothetical. 
However,  a  recent  NASA  study  suggests  that  space-borne  lidars  operating 
near  the  VUV  region  are  practical.  If  lidars  operating  Jn^  the  VUV  are 
developed  with  similar  capabilities  to  those  envisioned  by  NASA,  then 
significant  0*  measurements  can  be  made.  With  sufficient  output  power, 
lidars  have  an  extremely  good  inherent  range  resolution  of  about  5  m,  and 
horizontal  resolution  of  about  60  m.  Less  output  power  or  greater  range 
gives  less  resolution.  (The  NASA  instrument  has  a  range  resolution  of  1 
km  and  a  horizontal  resolution  of  10  km  in  a  typical  case.)  To  observe 
other  ionospheric  layers  would  require  a  lidar  capable  of  detecting  such 
molecular  ion  species  as  02*  and  NO'*'.  This  would  probably  require  a  second 
or  third  laser  transmitter  but  would  utilize  the  same  telescope  receiver. 

Thomson  scatter  radars  have  substantial  capabilities  for  observ¬ 
ing  the  ionosphere  from  the  ground  but  would  be  cumbersome  in  space. 
Because  the  modified  Thomson  scattering  cross  section  is  small,  both  a 
high  power  transmitter  and  a  large  antenna  are  required.  While  it  is 
within  current  technology  to  achieve  both  the  high  power  and  large  size  on 
spacecraft  such  as  the  Space  Shuttle,  a  space-borne  Thomson  scatter  radar 
would  still  be  limited  by  ground  and  auroral  clutter  problems  and  not  have 
a  particularly  good  resolution. 

Faraday  and  phase  methods  have  been  used  on  spacecraft  in  vari¬ 
ous  ways  for  two  decades.  These  methods  normally  require  a  cooperating 
ground  station  or  second  satellite  and  thus  lack  flexibility.  The  methods 
generally  yield  the  column  electron  density,  although  it  is  possible  to 
extract  the  peak  density  of  the  F-region  under  some  circumstances.  Faraday 
and  phase  methods  would  be  useful  as  part  of  a  larger  ionospheric  observa¬ 
tory,  but  probably  not  as  a  stand-alone  technique. 
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Many  passive  optical  or  chemiluminescent  methods  are  possible 
with  current  technology  and  relatively  simple  sensors.  For  instance,  a 
downward  looking  sensor  sampling  emissions  at  the  1356  A  and  6300  A  lines 
of  atomic  oxygen  can  be  used  to  provide  the  height  and  electron  density  of 
the  F-region  peak.  To  convert  the  column  emission  rates  that  the  sens'’’" 
measures  to  electron  density  requires  a  knowledge  of  chemistry  rate  coef¬ 
ficients,  neutral  densities,  and  the  shape  of  the  electron  density  pro¬ 
file.  A  limb  scanning  method  looking  for  the  same  emissions  can  provide 
the  shape  of  the  profile  as  well  as  the  electron  density  along  a  strip 
several  hundred  kilometers  wide  centered  on  the  orbital  plane  of  the 
satellite.  The  limb  scanning  photometer  measures  the  column  emission  rate 
but,  by  a  series  of  overlapping  measurements,  it  is  possible  to  extract 
the  volume  emission  rate  from  which  the  electron  density  can  be  calcu¬ 
lated.  Use  of  this  technique  in  conjunction  with  a  topside  sounder  reduces 
some  of  the  uncertainty  inherent  in  the  calculations. 

In  situ  methods  will  clearly  provide  the  most  detailed  measure¬ 
ments  of  the  ionosphere.  However,  by  definition,  they  have  little  areal 
survey  capability  without  a  great  many  spacecraft.  Therefore,  such  tech¬ 
niques  have  been  essentially  ignored  here. 

Bremsstrahlung  is  an  indirect  method  for  determining  ionospheric 
electron  density  which  is  potentially  useful  in  high-latitude  regions 
only.  By  measuring  the  bremsstrahlung  x-ray  flux  associated  with  a  dif¬ 
fuse  aurora,  it  may  be  possible  to  determine  the  flux  of  precipitating 
electrons.  This  flux  can,  in  turn,  be  related  to  the  ionospheric  electron 
density  in  regions  where  the  precipitating  flux  provides  the  dominant 
ionization.  Because  various  spectra  of  precipitating  electrons  will  pro¬ 
duce  the  same  bremsstrahlung,  it  is  necessary  to  select  a  most  probable 
spectra  based  on  a  given  set  of  bremsstrahlung  data.  An  approach  based  on 
a  maximum  entropy  formalism  seems  the  most  tractable.  The  lack  of  a 
unique  spectrun  may  not  cause  overly  great  difficulty  as  long  as  most  of 
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the  possible  spectra  produce  nearly  the  same  ionization.  This  work  is 
still  in  progress. 

Hybrid  methods,  incorporating  two  or  more  of  the  above  tech¬ 
niques,  would  clearly  be  superior  to  any  of  the  methods  alone  because  of 
the  complimentary  nature  of  the  various  techniques.  In  several  cases,  one 
method,  such  as  topside  sounding,  provides  the  calibrations  necessary  for 
another  method,  such  as  a  scanning  photometer  looking  for  atomic  oxygen. 
The  scanning  photometer  then  provides  the  electron  density  in  the  region 
adjacent  to  the  subsatellite  track,  which  cannot  be  seen  by  the  sounder. 
Similarly,  a  laser  radar  would  benefit  from  cross  calibration  with  a  top¬ 
side  sounder  and  would  then  be  able  to  accurately  survey  a  far  larger  area 
than  could  the  sounder.  The  same  argument  applies  to  a  topside  sounder 
working  with  a  Thomson  scatter  radar,  a  Faraday  and  phase  technique,  or  a 
bremsstrahlung  method. 

A  space-borne  ionospheric  observatory  consisting  of  a  topside 
sounder  such  as  the  Digisonde  and  a  laser  radar  would  have  substantial 
capability  for  performing  an  areal  survey,  a  capability  far  beyond  any 
currently  available.  For  little  more  effort,  complexity,  or  cost,  one 
could  include  a  scanning  photometer  and  bremsstrahlung  x-ray  detector. 
This  would  provide  considerable  redundancy  enabling  substantial  cross 
calibration.  Not  only  would  such  an  instrument  provide  an  excellent  areal 
survey  of  ionospheric  electron  density  over  most  of  the  globe  day  or 
night,  but  also  additional  information  useful  in  ionospheric  studies. 

Specific  recommendations  for  the  future  include: 

1)  Continued  pursuit  of  laser  radars. 

2)  Continued  pursuit  of  bremsstrahlung  techniques. 

3)  Design  of  a  hybrid  spacecraft. 
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SECTION  2 
INTRODUCTION 


Although  the  earth's  ionosphere  has  been  studied  in  many  differ¬ 
ent  ways  for  over  fifty  years,  it  is  still  pursued  for  scientific  as  well 
as  operational  reasons.  The  gross  structure  of  the  ionosphere  is  well 
understood,  while  the  detailed  structure  and  the  evolution  of  the  struc¬ 
ture  with  time  are  only  partially  understood.  In  addition,  the  incoming 
electromagnetic  and  particulate  radiations  which  influence  and  maintain 
the  ionosphere  are  variable  and  incompletely  known.  As  a  consequence, 
predictions  of  ionospheric  behavior  for  such  applications  as  HF  radio 
propagation  are  possible  but  subject  to  error.  Hence,  the  ionosphere 
remains  a  subject  of  substantial  interest. 

This  paper  concentrates  on  techniques  for  directly  or  indirectly 
measuring  ionospheric  electron  density  from  a  satellite  platform.  Direct 
measurements  of  electron  density  are  important  to  those  concerned  with 
radio  propagation  as  well  as  to  those  building  or  verifying  models. 
Indirect  measurements  of  electron  density  such  as  measurements  of  incoming 
radiation  can  be  used  to  calculate  the  electron  density  using  large  com¬ 
puter  codes  which  incorporate  the  basic  physics  of  the  ionosphere. 

Satellite  measurements  are  emphasized  by  this  study  so  that 
large  or  remote  regions  of  the  ionosphere  can  be  observed.  Because  cur¬ 
rently  existing  satellites  can  only  measure  electron  density  over  a  very 
limited  region  of  the  ionosphere  during  any  given  orbit,  this  investiga¬ 
tion  looks  for  methods  by  which  a  more  general  areal  survey  can  be  made. 
Satellite  instruments  have  a  large  areal  coverage  but  are  necessarily 
limited  by  space,  weight,  power,  and  repair  considerations  to  less 
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elaborate  or  tempermental  instruments.  Most  of  the  techniques  discussed 
here  involve  self-contained  satellite  instruments  but  a  few  require  a 
cooperating  ground  station. 


Of  the  many  ionospheric  parameters  which  can  be  and  have  been 
measured,  electron  density  is  the  most  fundamental  in  scientific  studies 
and  the  most  important  to  those  attempting  to  use  the  ionosphere  for  com¬ 
munications.  Other  ionospheric  parameters  such  as  positive  and  negative 
ion  species  and  concentrations,  neutral  species  and  concentrations,  and 
ion  and  electron  temperatures  are  important  but  not  of  primary  interest 
here. 


Because  it  is  difficult  to  measure  all  aspects  of  ionospheric 
electron  density  simultaneously,  particular  aspects  of  this  density  are 
emphasized  at  a  given  time  depending  on  the  particular  application.  For 
instance,  electron  density  is  pertinent  to  HF  propagation,  while  integrat¬ 
ed  electron  density  (column  electron  density  or  total  electron  content)  is 
useful  to  satellite  communications.  Similarly,  small  scale  electron 
density  structure  is  interesting  in  some  contexts,  while  large  scale 
structure  is  imporcant  in  others.  Sometimes,  three  dimensional  spatial 
variation  of  electron  density  is  wanted,  while  under  other  circumstances 
time  development  in  one  region  is  important.  If  one  technique  could  be 
developed  to  measure  all  of  these  aspects  at  once,  the  problem  would  be 
solved.  Unfortunately,  no  such  technique  exists  or  is  likely  to  exist. 
Some  techniques  can  measure  more  than  one  aspect  of  electron  density  but 
are  not  good  for  others.  For  instance,  the  ISIS  satellites  have  the  best 
existing  capability  for  measuring  ionospheric  electron  density.  The 
topside  sounders  on  these  satellites  measure  electron  density  along  the 
sub-satellite  track.  They  obtain  the  electron  density  in  two  dimensions 
but  not  three.  They  see  only  the  topside  and  not  the  bottoms ide.  They 
measure  the  spatial  variation  of  electron  density  in  a  north-south  plane 
but  not  an  east-west  plane.  They  do  not  remain  long  enough  in  any  region 
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to  watch  the  time  evolution  of  the  electron  density.  Despite  these  short¬ 
comings,  the  ISIS  topside  sounders  have  collected  vast  quantities  of  high¬ 
ly  useful  data.  Some  of  the  shortcomings  of  various  instruments  can  be 
overcome  by  using  two  or  more  instruments  with  complimentary  capabilities. 
Although  a  technique  without  shortcomings  would  be  preferred,  it  is  pro¬ 
ductive  to  consider  any  technique  which  may  yield  a  different  aspect  of 
electron  density. 

This  paper  considers  a  broad  range  of  techniques  from  those  that 
have  been  tried  and  proven  to  those  that  are  more  speculative.  Old  tech¬ 
niques  have  been  reexamined  to  look  for  additional  potential,  while  sever¬ 
al  new  techniques  have  been  investigated  to  the  point  where  their  current 
potential  can  be  evaluated.  Active  methods  which  probe  the  ionosphere 
with  radio  waves  or  light  are  discussed  first,  while  passive  methods  which 
use  emissions  from  the  ionosphere  are  discussed  last.  Several  of  the 
active  techniques  directly  yield  electron  density  while  most  of  the  pas¬ 
sive  techniques  yield  the  less  useful  column  electron  density. 

Heading  the  list  of  satellite  techniques  is  topside  sounding, 
because  this  technique  is  the  best  one  in  current  use  for  measuring  elec¬ 
tron  density  over  wide  areas.  It  is  the  technique  against  which  others 
must  be  measured.  Although  topside  sounding  has  been  used  for  nearly 
twenty  years,  it  has  additional  potential  which  is  discussed.  Two  new 
techniques,  laser  radars  and  Thompson  scatter  radars,  offer  the  possibil¬ 
ity  of  improved  area  coverage  and  improved  electron  density  measurements. 
Various  Faraday  and  phase  techniques  are  discussed.  These  have  limited 
potential  because  they  require  a  cooperating  ground  station  or  second 
satellite  and  still  yield  only  column  electron  density.  Many  chemilumi¬ 
nescent  methods  are  possible,  some  of  which  have  considerable  potential. 
In  situ  measurements  have  the  virtue  of  high  resolution  but  the  drawback 
of  very  limited  areal  coverage.  Ionospheric  electron  density  in  the  polar 
regions  is  indirectly  related  to  bremsstrahlung  x-rays  from  precipitating 
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electrons.  Measurements  of  this  radiation  can  provide  information  on  an 
important  source  of  ionization  in  the  polar  ionosphere.  To  derive  elec¬ 
tron  density  from  such  measurements  requires  considerable  calculation. 
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SECTION  3 
TOPSIDE  SOUNDERS 


The  satellite  sensing  technique  currently  offering  by  far  the 
best  areal  survey  of  ionospheric  electron  density  is  topside  sounding.^ 
The  two  ISIS  satellites^  built  by  the  Communications  Research  Centre  in 
Ottawa,  Ontario,  Canada  have  topside  sounders  that  have  been  operating  in 
polar  orbits  for  more  than  ten  years.  Although  these  instruments  do  an 
excellent  job  of  sounding,  advances  in  technology  over  the  last  decade  now 
permit  considerably  more  sophisticated  measurements.  An  instrument  called 

the  "Digisonde"^ built  by  the  University  of  Lowell  in  Lowell,  Massa¬ 
chusetts,  incorporates  the  new  technology  and  could  rather  easily  be  modi¬ 
fied  for  spacecraft  applications. 

Topside  sounding  in  its  simplest  form  is  analogous  to  the  sonar 
technique  used  by  ships  to  probe  the  ocean.  The  topside  sounder  broad¬ 
casts  a  short  radio  wave  pulse  at  a  frequency  anywhere  from  about  .1  MHz 
up  to  10  or  20  MHz  and  listens  for  a  return.  The  pulse  reflects  off  of 
the  ionosphere  below  the  satellite  when  the  pulse  encounters  an  electron 
density  with  a  plasma  frequency  equal  to  the  radio  wave  frequency.  This 
electron  density  is  given  by  the  well  known  relationship 

N  =  1.24  X  10**  f2  (3-1) 

P 

where  N  is  the  electron  density  in  electrons/cm^  and  fp  is  the  plasma 
frequency  in  Wz.  By  measuring  the  time  it  takes  for  the  pulse  to  return 
to  the  satellite,  the  virtual  distance  from  the  satellite  to  a  particular 
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electron  density  in  the  ionosphere  can  be  found.  If  a  succession  of  radio 
wave  pulses  are  used  starting  near  0.1  MHz  and  proceeding  in  steps  up  to 
10  or  20  MHz,  the  vertical  profile  of  the  topside  ionosphere  below  the 
satellite  can  be  obtained. 

The  satellite  does  not,  however,  see  all  of  the  ionosphere 
beneath  it.  The  bottoms ide  cannot  be  seen  with  this  technique  because 
pulses  with  a  sufficiently  high  frequency  to  penetrate  through  the  topside 
are  not  reflected  by  the  lower  electron  densities  existing  below  the  peak 
of  the  F-layer.  Hence,  the  bottomside  F2,  the  FI,  E,  and  D  layers  go 
undetected  by  this  method.  Also,  the  topside  sounder  normally  receives 
the  reflected  echo  from  directly  beneath  the  satellite  even  though  the 
sounder  transmitting  antenna  is  essentially  isotropic.  This  occurs  because 
the  smooth  ionosphere  reflects  radio  waves  as  a  mirror  reflects  light. 
Because  of  the  curvature  of  the  earth  and  the  corresponding  curvature  of 
the  ionosphere,  only  the  radio  wave  reflection  occurring  beneath  the 
satellite  returns  a  signal  to  the  satellite.  This  means  that  topside 
sounders  can,  in  general,  only  scan  the  topside  ionosphere  along  the  sub¬ 
satellite  track.  As  the  satellite  moves  in  its  orbit  around  the  earth,  it 
obtains  an  excellent  two-dimensional  view  of  the  topside  but  misses  three 
dimensional  characteristics,  misses  the  bottomside,  and  misses  any  varia¬ 
tion  with  time.  Time  variations  and  three-dimensional  characteristics 
can,  to  some  extent,  be  pieced  together  using  statistical  studies  of  many 
passes. 

The  two  ISIS  spacecraft  accomplish  topside  sounding  and  eight  to 
ten  other  measurements  of  the  ionospheric  parameters.  They  are  moderate 
size,  moderate  weight  instruments,  consume  only  125  watts  in  full  opera¬ 
tion,  and  are  highly  reliable.  After  more  than  ten  years  of  operation, 
most  of  the  major  spacecraft  systems  in  both  satellites  are  still  func¬ 
tioning  (including  the  sounders).  Table  2  summarizes  the  characteristics 
of  ISIS  II.  These  characteristics  are  similar  to  those  of  ISIS  I. 


Table  2.  Sunmary  of  ISIS  II  spacecraft  characteristics  (adapted 
from  Reference  2). 

SPACECRAFT  SIZE 
Shape 
Height 
Diameter 
Weight 

ORBIT 

Approximately  1400  km  circular 
Inclination  88.7*  prograde 

STABILIZATION 

Spin  stabilized 

ALTITUDE  SENSING 

Six-probe  flux-gate  magnetometer 
Digital  solar  aspect  sensor 

SPIN  AND  ATTITUDE  CONTROL  SYSTEM 

Normal  spin  rate  is  near  3  rpm  with  a  capability  for  changing 
the  spin  axis  into  either  of  two  configurations 

a.  Spin  axis  in  the  orbit  plane 

b.  Spin  axis  in  cartwheel  configuration 

POWER  SYSTEM 

11,008  solar  cells 
3  nickel -cadmium  batteries 
5  main  dc  to  dc  converters 

125  watts  for  4  hours/day  operation  in  70  percent  sunlight 
2  consecutive  pole-to-pole  passes 


approximates  an  oblate  spheroid 
1.19  m 
1.27  m 
261  kg 


Table  2  (continued).  Summary  of  ISIS  II  spacecraft  characteristics 

(adapted  from  Reference  2). 


COMMAND  SYSTEM 

Multiple  digital  tone  AM  system 
2  command  receivers  connected  redundantly 
216  command  capability 


TELEMETRY  SYSTEM 

Transmitter  #1 

Frequency  136.080  MHz  with  100  kHz  bandwidth 

Power  4  watts 

Modulation  FM/PAM/PDM  (housekeeping  data) 

FM/PCM  (Auroral  photometer) 


Transmitter  #2 

Frequency  136,59  MHz  with  50  kHz  bandwidth 

Power  2  watts 

Modulation  PM/PCM  (experiments  and  housekeeping) 


Transmitter  #3 

Frequency  401.75  MHz  with  300  kHz  bandwidth 

Power  4  watts 

Modulation  FM  and  PCM 
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The  sounders  on  board  these  spacecraft  can  be  operated  in  either 
a  swept  frequency  or  fixed  frequency  mode.  They  have  a  transmitter  pulse 
power  of  400  watts  at  a  pulse  repetition  frequency  of  45  pps  and  a  pulse 
width  of  100  liSec.  This  pulse  width  determines  the  inherent  vertical 
resolution  of  15  km.  The  inherent  horizontal  resolution  is  determined  by 
the  size  of  the  first  Fresnel  zone  of  the  reflecting  surface  beneath  t^p 
satellite.  This  is  typically  5  km  but  depends  on  the  range  of  reflecting 
surface  from  the  satellite  and  the  sounding  frequency.  Irregularities 
with  a  smaller  scale  size  than  5  km  are  observable  under  some  special 
circumstances.  One  example  is  given  below.  The  sounder  transmitter  feeds 
the  sounder  ante,''pas  via  a  passive  network  which  provides  cross-over 
coupling  as  a  function  of  frequency  between  the  long  (73  m)  and  short  (19 
m)  dipole  anternas.  T,ie  analog  receiver  uses  the  same  antenna  system  and 
operates  on  one  ’  -.everal  different  broad  bands.  Table  3  summarizes  the 
characteristics  of  the  ISIS  II  sounder  while  Figure  1  provides  a  block 
diagram  of  the  sounder  system.  For  completeness.  Table  4  is  included  to 
summarize  ^he  otner  ionospheric  sensors  on  ISIS  II.  Again,  ISIS  I  is  very 
similar  to  ISIS  II. 

In  reality,  topside  sounding  is  considerably  more  complex  than 
it  may  seem.  Figure  2  shows  perhaps  the  simplest  of  the  topside  ionograms 
obtained  by  ISIS.  This  particular  ionogram  comes  from  ISIS  II  and  was 
obtained  in  1978  on  Kwajalein  Atoll. ^  It  shows  various  sounder  traces 
plotted  as  a  function  of  sounding  frequency  and  apparent  range  from  the 
satellite.  The  apparent  range  has  been  computed  from  the  measured  time 
delay  by  multiplying  by  the  speed  of  light  and  dividing  by  two.  The  width 
of  the  X  and  0  traces  show  the  15  km  inherent  vertical  resolution  defined 
by  the  transmitted  pulse  width.  Separate  X,  0,  and  Z  traces  occur  because 
the  ionosphere  is  imbedded  in  the  earth's  magnetic  field.  All  of  these 
traces  can  be  explained  by  the  Appleton-Hartree  formalism.^  The  various 
resonance  spikes  occur  because  the  satellite  is  within  the  ionosphere. 


Table  3.  ISIS  II  topside  sounder  characteristics  (adapted  fron 
Reference  2). 


SWEPT  FREQUENCY  SOUNDER 
Power  output 
Pulse  width 
Antennas 


400  watts  at  PRF  45  pps 
100  usec 

Two  crossed  dipoles  73  tn  and  19  m  at 
tip  to  tip 


Frequency  range 

Normal  sweep  to  0.1  to  10  Wiz  in  10  sec. 
Extended  sweep  0.1  to  20  MHz  in  20  sec. 
Frequency  range  Sweep  rate 

0.1  to  2.0  MHz  0.375  MHz/sec. 

2.0  to  5.0  MHz  1.125  MHz/sec. 

5.0  to  20.0  MHz  1.5  MHz/sec. 


Frequency  markers 

0.1,  0.25,  0.5,  0.75,  1.0,  1.25,  1.5,  1.75,  2.0,  3.0,  4.0, 
5.0,  6.0,  7.0,  8.0,  9.0,  10.0,  12.0,  14.0,  16.0,  18.0,  20.0 

Frequency  resolution  between  markers  is  ±  8  kHz. 

Range  marker 

11.11  msec  after  zero  range  mark  (not  on  ISIS  I) 

Sounder  receiver  bands 

(a)  0.1  -  1.0  MHz 

(b)  1.0  -  2.0  MHz 

(c)  2,0  -  3.0  MHz 

(d)  3.0  -  5.0  Mhz 

(e)  5.0  -  8.0  MHz 

(f)  8.0  -  13.0  MHz 

(g)  13.0  -  20.0  MHz 


Table  3  (continued).  ISIS  II  topside  sounder  characteristics 

(adapted  from  Reference  2). 

FIXED  FREQUENCY  SOUNDER 

6  fixed  frequencies  are  available  by  command: 

(0.12,  0.48,  1.0,  1.95,  4.0,  and  9.303  MHz) 

These  occur  either  during  the  flyback  period  of  the  sounder 
sweep  or  for  a  full  sounder  frame  period  on  alternate 
frames. 


Block  diagram  of  Sounder  System  (from  Reference  2) 


Table  4.  Sunmiary  of  other  experiments  on  ISIS  II  (adapted  from 
Reference  2). 

VLF  EXPERIMENT 

Receiver  for  radio  noise  studies  (.05  to  30  Khz) 

Optical  onboard  receiver  for  ion  resonance  studies 

ENERGETIC  PARTICLE  DETECTORS 

Detectors  to  measure  electron  and  proton  intensity,  angular 
distributions,  and  energy  spectra 

SOFT  PARTICLE  SPECTROMETER 

Electrostatic  analyzer  to  measure  directional  intensity  and 
energy  spectra  of  protons  and  electrons  from  10  eV  to  10  keV 

ION  MASS  SPECTROMETER 

Magnetic  spectrometer  to  measure  positive  ion  mass  (1-64  amu) 
LANGMUIR  PROBE 

Two  cylindrical  electrostatic  probes  at  opposite  ends  of  the 
spacecraft  to  mesure  electron  density  and  temperature 

RETARDING  POTENTIAL  ANALYZER 

Direct  measurement  of  electron  temperature,  ion  temperature,  ion 
composition,  and  charged  particle  density 

BEACON 

Two  100  mW  beacons  to  investigate  ionospheric  irregularities 
(136.410  and  137.950  MHz) 


Table  4  (continued).  Summary  of  other  experiments  on  ISIS  II 

(adapted  from  Reference  2). 

COSMIC  NOISE 

Monitors  background  radio  noise  levels  due  to  galactic,  solar, 
and  ionospheric  sources  using  A6C  of  sounder  receiver  (0.1  to  16 
MHz) 

ATOMIC  OXYGEN  RED  LINE  PHOTOMETER  (NOT  ON  ISIS  I) 

Maps  global  distribution  in  the  intensity  of  the  6300  A  emission 
from  atomic  energy 

AURORAL  SCANNER  (NOT  ON  ISIS  I) 

Maps  the  distribution  of  auroral  emissions  at  5577  and  3914  A  on 
the  dark  side  of  the  earth 


Resonance 


daytime  ionogram®. 


tven  though  some  of  the  resonance  s()ikes  extend  to  a  large  apparent  range, 
they  only  represent  resonances  occurring  in  the  immediate  vicinity  of  the 
satellite.  Similarly,  the  apparent  ranges  of  the  X,  0,  and  Z  traces  are 
not  equal  to  their  real  ranges,  although  the  traces  do  represent 
reflections  from  below  the  satellite.  For  these  traces,  the  apparent 
range  is  always  greater  than  the  real  range  because  the  sounder  pulse 
always  sustains  some  retardation  in  phase  velocity  just  before 
reflection.  Such  retardation  is  well  understood,  so  that  topside 
ionograms  can  be  rather  easily  reduced  to  true  height  profiles.® 

Not  all,  or  even  most,  ionograms  are  as  simple  or  as  readily 
understood  as  the  one  in  Figure  2.  For  instance.  Figure  3  shows  an  iono- 
grani  with  a  succession  of  traces  at  progressively  greater  ranges.  In 
keeping  with  their  appearance,  they  are  labeled  "Multiple  Echoes."  The 
strong  multiple  echo  at  greatest  range  is  called  a  "Conjugate  Echo"  for 
reasons  explained  below.  The  departure  of  Figure  3  from  the  simple  form 
of  Figure  1  occurs  because  the  ionosphere  is  not  smooth.  Returns  in 
Figure  2  are  coming  from  both  the  sub-satellite  track  as  well  as  from 
other  locations  relatively  far  removed  from  directly  beneath  the  satel¬ 
lite.  In  this  case,  the  satellite  has  a  far  reaching  view.  It  turns  out 
that  these  multiple  echoes  are  ducted  echoes.®*^®  They  are  evidence  that 
magnetic  field  aligned  tubes  of  slightly  depleted  electron  density,  per¬ 
haps  1  km  in  diameter,  exist  below  the  satellite  and  extend  for  hundreds 
or  thousands  of  kilometers  along  the  field.  They  differ  from  the  sur¬ 
rounding  electron  density  by  a  few  percent  and  have  widths  of  a  few  kilo¬ 
meters.  Tiie  conjugate  echo  is  a  special  ducted  echo  where  the  satellite 
is  actually  inside  the  duct.  Radio  waves  trapped  within  this  duct  can 
echo  back  and  forth  from  one  magnetic  conjugate  region  to  the  other  many 
times  giving  a  strong  return  each  time  they  pass  the  satellite. 

The  ionogram  shown  in  Figure  4  illustrates  another  basic  type  of 
disturbed  ionosphere  seen  by  ISIS.  This  is  one  version  of  the  classic 
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"spread-F”  often  seen  on  ionogranis.  It  is  evidence  of  a  great  deal  of 
small  scale  structure  below  the  satellite.  Spread-F  comes  in  several 
forms.  This  particular  one,  known  as  an  "exploded  trace,"  occurs  only  in 
the  equatorial  ionosphere  and,  in  fact,  only  within  the  equatorial  anomaly 
region  from  lb°N  to  Ib^S  magnetic  latitude.  Recently,  such  traces  have 
been  correlated  with  radar  backscatter  plumes  and  severe  satellite  signal 
scintillations.^  They  are  probably  caused  by  the  equatorial  bubbles  under 
intensive  investigation  for  the  last  few  years. 

As  with  many  complex  topside  ionograms,  some  controversy  exists 
as  to  whether  the  spread  returns  in  Figure  4  represent  overdense  or  under- 
dense  reflections.  Reference  12  views  exploded  traces  as  an  overdense 
layer  of  irregularities  above  the  otherwise  smooth  ionosphere.  This  layer 
has  a  limited  east-west  but  long  north-south  extent.  The  irregularities 
supposedly  produce  overdense  reflections  but  are  sufficiently  sparse  that 
some  of  the  sounder  signal  can  penetrate  to  the  ionosphere  below.  In 
contrast.  Reference  13  views  the  irregularities  as  underdense,  producing 
partial  radio  wave  reflections  because  their  cross  sections  are  small 
compared  to  the  wavelengths  of  the  sounder  radio  waves.  The  irregularities 
then  need  not  differ  greatly  in  density  from  the  surrounding  plasma.  They 
can  be  either  slight  enhancements  or  depletions.  For  various  reasons, 
underdense  reflection  appears  the  more  logical  choice.® 

Although  topside  ionograms  normally  provide  only  a  two-dimen¬ 
sional  view  of  the  topside  ionosphere,  sometimes  a  three-dimensional  view 
is  obtained.  Figure  5  shows  two  exploded  traces.  One  trace  is  attached 
to  the  smooth  primary  trace  and  thus  represents  a  plume  beneath  the  satel¬ 
lite  while  the  other  represents  a  separate  plume  displaced  to  the  east  or 
west  of  the  satellite.  From  successive  ionograms,  it  is  clear  that  these 
plumes  are  indeed  separated  east-west  and  have  a  very  long  north-south 
extent.  Aspect  sensitive  underdense  reflection  of  the  radio  waves  from 
the  irregularities  permits  the  east-west  extent  of  the  irregularities  to 
be  seen  in  a  given  ionogram  while  suppressing  the  north-south  extent. 
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ISIS  II  1977  Day  228  1021  Z  2109  LT  (6’N.  162*E.  1415  km  high). 
Exploded  trace  and  remote  spread®. 


As  the  satellite  travels  to  the  north,  successive  ionograms  provide  the 
three  dimensional  picture. 

This  three-dimensional  picture  has  definite  ambiguities.  For 
instance,  there  is  no  way  to  tell  whether  the  remote  echo  is  coining  from 
the  east  or  from  the  west.  Also,  it  would  be  interesting  to  know  if  the 
irregularities  are  stationary  or  are  rising  as  equatorial  bubbles  are 
thought  to  rise.  The  ISIS  topside  sounders  are  unable  to  answer  these 
questions  even  though  the  information  could  be  obtained  from  the  reflected 
sounder  pulses. 

A  new  instrument  called  the  "Digi sonde"  developed  over  the  last 
decade  at  the  University  of  Lowell  attempts  to  use  all  of  the  information 
available  from  the  reflected  wave  and  thereby  eliminate  some  ambiguities 
that  have  previously  existed  with  sounding  techniques.  The  parameters 
that  can  be  determined  from  the  reflected  wave  are: 

1.  Range  (group  travel  time) 

2.  Amplitude 

3.  Phase 

4.  Doppler  offset  from  transmitted  frequency 

5.  Angle  of  arrival 

6.  Wave  polarization  (0  or  X) 

7.  Curvature  of  wave  front. 

The  ISIS  sounders  can  only  measure  the  first  parameter  in  a  quantitative 
fashion,  while  the  Digisonde  can  measure  the  first  six  quantitatively. 
The  last  three  parameters  require  an  antenna  array  which  is  more  practical 
for  a  bottomside  sounder  than  for  a  topside  sounder. 

Digital  techniques  are  the  key  to  the  Digi sonde's  success.  Such 
techniques  have  become  feasible  in  the  last  decade  with  the  advent  of 
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microprocessors  and  associated  digital  integrated  circuits.  The  techniques 
used  by  the  Digisonde  to  measure  the  parameters  listed  are  summarized  in 
Table  5. 

Although  it  is  possible  to  monitor  simultaneously  the  entire 
five-dimensional  complex  amplitude  space  (range,  frequency,  Doppler, 
azimuth,  and  zenith  angle),  this  is  neither  feasible  nor  desirable.  The 
quantity  of  raw  data  would  be  excessive.  On  the  other  hand,  it  is  impor¬ 
tant  to  record  as  much  raw  data  as  possible  and  not  immediately  process  it 
into  desired  results.  The  Digisonde  avoids  this  dilemma  by  digitally 
preprocessing  the  data  to  reduce  redundancies  and  to  flag  each  data  point 
with  a  parameter  identification.  It  also  operates  in  one  of  two  modes: 
ionogram  or  drift  mode.  These  modes  are  complimentary  because  together 
they  contain  all  the  relevant  information  of  the  five-dimensional  space. 
In  the  ionogram  mode,  all  frequencies  and  range  bins  are  sampled  with  full 
amplitude  resolution  but  with  limited  information  on  incidence  angles  and 
Doppler  frequencies.  The  drift  mode  examines  only  a  few  frequencies  and 
range  bins  but  obtains  full  amplitude,  incidence  angle,  and  Doppler  infor¬ 
mation.  In  practice,  the  ionogram  mode  is  only  used  occasionally  to 
establish  the  best  frequencies  and  range  bins  for  the  Doppler  mode.  The 
Doppler  mode  then  provides  all  of  the  data  that  is  normally  required.  The 
multiplexed  output  data  still  retains  its  "raw"  character  but  has  been 
reduced  to  manageable  proportions. 

The  Digisonde  is  able  to  reduce  significantly  the  effects  of 
noise  interference  by  phase  coding  of  the  transmitted  signal.  When  used 
as  a  bottomside  sounder,  this  feature  greatly  improves  the  quality  of  the 
data  which  would  otherwise  be  poor  due  to  excessive  man-made  noise.  On 
the  topside  of  the  ionosphere,  man-made  noise  is  much  less  of  a  problem 
because  the  bottomside  reflects  most  of  it.  However,  on  some  occasions, 
especially  with  a  highly  disturbed  ionosphere,  man-made  noise  can  be  a 
problem  on  the  topside  also.  Under  these  circumstances,  the  Digisonde 
would  perform  better  than  ISIS. 


Table  5.  Oigisonde  measurement  techniques  (adapted  from  Reference  3). 


WAVE  MEASUREMENT 

PARAMETER  TECHNIQUE 


(1)  Range 


Timing  of  interval  between 
transmitted  and  received  pulse 


(2)  Amplitude 


(3)  Phase 


(4)  Doppler 


Coherent  integration  of  16  to 
256  quadrature  samples,  log 
compression  of  signal,  and 
digitization  in  1/2  dB  steps 

Ratio  of  integrated  quadrature 
S6mples.  Synchronization  of 
IF  and  sampling  time 

Discrete  complex  Fourier 
transform 


(5)  Angle  of  Arrival  Receiving  antenna  array  and 

antenna  switch 


(6)  Wave  Polarization  Polarized  receiving  antennas 

and  antenna  switch 
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The  Digisonde  is  a  versatile  instrument  with  a  wide  range  of 
programmable  features.  Table  6  details  these  while  Figure  6  provides  a 
block  diagram  of  the  equipment  as  it  is  configured  for  bottoms ide  sound¬ 
ing. 


Table  6  shows  that  the  basic  vertical  and  horizontal  resolution 
of  the  Digisonde  should  be  similar  to  that  of  the  ISIS  sounders.  Indeed 
the  pulse  widths  of  50,  100,  or  150  psec  are  close  to  those  of  ISIS  imply¬ 
ing  an  inherent  resolution  of  7.5  to  22.5  km.  However,  the  use  of  phase 
measurements  increases  the  group  height  resolution  significantly.  This 
group  height  is  given  by 

h'  =  ^  (3-2) 

Att  Af 

where  c  is  the  speed  of  light  and  A4>  is  the  change  in  phase  with  a  change 
in  frequency,  Af.  Item  number  3  in  Table  6  shows  that  a  dual  step  size 
for  frequency  can  be  selected.  This  is  for  use  with  phase  measurements. 
The  large  frequency  step  provides  high  resolution  while  the  small  step 
resolves  2-n  ambiguities  in  the  phase.  For  the  case  of  10  and  90  kHz  dual 
frequency  increments,  the  resolution  in  virtual  height  would  be 

Ah'  «  100  m  (3-3) 

which  is  two  orders  of  magnitude  better  than  ISIS. 

The  horizontal  resolution  and  field  of  view  of  the  Digisonde 
remain  similar  to  ISIS.  However,  the  Digisonde  transmits  a  10  kilowatt 
pulse  in  comparison  to  400  watts  for  ISIS.  This  14  dB  greater  power 
output  is  highly  desirable,  because  ISIS  often  looses  large  portions  of 
the  topside  X  and  0  traces  due  to  inadequate  power.  The  continuous  power 
required  by  the  Digisonde  is  typically  200  watts,  which  is  certainly 
within  the  realm  of  current  spacecraft  capabilities. 
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Table  6.  Progranmable  ionogram  parameters  (adapted  from  Reference  4) 


■ 

Function 

Dimension 

Va  1  ue 

1 

Beg  In  frequency 

HMz 

0,  1,  2,  ....  39 

2 

End  frequency 

MHz 

1  p  2p  •  •  •  f  40 

3 

Frequency  Increments 

kHz 

5,  10,  25.  50,  100,  200 

5  &  20,  10  &  40,  10  &  90.  20  &  180 

4 

Range  Increments 

km 

0.5,  1.0,  1.5,  2.0,  3.0,  4.5,  5.0,  6.0,  10.5 

5 

Range  samples 

128,  256 

6 

Doppler  resolution 

Hz 

0.2,  0.4,  0.8,  1.6,  3.2,  6.3,  12.5,  25.0* 

7 

Number  of  Doppler  lines 

1,  2,  4,  8,  16,  64,  128 

8 

4mpl Itude  resolution 

dB 

0.5,  1.0 

9 

Phase 

deg 

22.5 

10 

Pulse  repetition  rate 

Hz 

50,  100,  200,  400 

11 

Pul se  width 

US 

50,  100,  150 

12 

SampI  es  per  frequency 

16,  32,  64,  128  ,  256 

13 

Azimuth  resolution 

deg 

30,  60,  90,  120,  180,  none 

14 

Wave  polarization 

ordinary,  extraordinary,  ordinary  and  extraordinary 

15 

Range  start 

km 

7.5,  37.5,  60.0,  82.5,  112.  5,  172.5,  232.5,  352.5,  832.5,  712.5 

16 

Range  window 

Rj  1  se  per  1  od 

-4,  -3,  -2,  -1,  no  Ident. ,  +1,  +2,  t3,  +4,  +5 

17 

Receiver  gain  adj. 

dB 

0,  6,  12,  18  ,  24  ,  30,  36,  42 

18 

Ionogram  rate 

hr-1 

1  to  120 

*Dfv1de  by  nunber  of  antenna  configurations  used. 
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Figure  6.  Digisonde  system  diagram  (from  Reference  4) 


In  terms  of  readily  available  technology  that  could  be  adapted 
to  spacecraft  use  and  would  provide  a  significant  improvement  in  our 
capability  to  monitor  ionospheric  electron  density,  the  Digisonde  is  a 
clear  first  choice.  It  is  our  understanding  that  plans  are  proceeding  to 
put  several  Digisonde-1  ike  instruments  into  polar  orbits  by  the  end  of 
this  decade. 
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SECTION  4 

LASER  RADAR  (LIDARS) 


Laser  radars  or  "lidars"  appear  to  have  the  greatest  long-term 
potential  for  high-quality  observations  of  the  ionosphere  from  a  satel¬ 
lite.  Some  lidar  observations  of  the  ionosphere  have  been  made  with 
ground  based  instruments but  these  have  been  limited  to  measurements 
of  neutral  metal  atoms  in  the  lower  ionosphere.  While  these  measurements 
are  interesting  in  their  own  right,  they  do  not  relate  directly  to 
ionospheric  electron  density.  It  is  necessary  to  measure  the 
concentrations  of  the  major  ionic  constituents  (especially  0'*')  to  obtain 
electron  density.  Such  measurements  require  a  tunable  lidar  operating  in 
the  vacuuti  ultraviolet  (VUV)  near  80  nm  (for  0'^)  instead  of  589  nm  (for 
Na) .  Currently,  no  tunable  lidar  exists  with  a  wavelength  below  100  nm, 
but  they  are  certainly  possible.  If  one  were  developed,  it  would  be 
ideally  suited  to  a  spacecraft  instead  of  a  ground-based  application, 
because  the  lower  atmosphere  is  opaque  in  the  VUV. 

A  laser  radar  operates  in  much  the  same  way  as  any  other  radar. 
It  emits  a  short  pulse  of  radiation  (in  this  case,  light)  and  waits  for 
the  pulse  to  be  scattered  back  by  the  mediun  being  probed.  By  measuring 
the  intensity  of  the  returning  signal  as  a  function  of  the  time  it  takes 
to  return,  it  is  possible  to  determine  the  density  of  the  scattering 
medium  as  a  function  of  distance  from  the  lidar. 

Figure  7  illustrates  in  more  detail  how  the  lidar  technique 
would  be  applied  from  a  spacecraft.  The  short  high-power  pulse  of  light 
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FRECEDIN9  PAOK  BUMCp^OV  FILMD) 


Illustration  of  a  Ildar  observing  the  ionosphere  (adapted  from 
Reference  22). 


from  the  laser  transmitter  is  directed  at  the  earth's  ionosphere.  This 
pulse  of  light  has  the  appropriate  wavelength  to  be  scattered  by  the  iono¬ 
spheric  constituent  of  interest.  A  telescope  coupled  to  a  photomultiplier 
tube  detector  receives  scattered  light  from  the  laser  as  well  as  scattered 
light  (noise)  from  other  sources  such  as  the  sun.  This  light  is  converted 
to  electrical  pulses  by  a  photomultiplier  tube  (or  photodiode)  and  subse¬ 
quently  processed  by  on-board  electronics. 

Lidars  offer  the  possibility  of  extremely  good  spatial  resolu¬ 
tion.  Because  the  typical  lidar  pulse  length  is  only  30  ns,  a  range 
resolution  of  less  than  5  m  is  possible.  This  compares  with  the  inherent 
range  resolution  of  15  km  for  the  ISIS  topside  sounders.  The  horizontal 
resolution  of  a  lidar  is  also  high.  The  typical  divergence  of  a  laser 
pulse  is  1  mrad  which  produces  a  spread  of  only  300  m  at  a  range  of  300 
km.  By  restricting  the  field  of  view  of  the  receiving  telescope  to  .2 
mrad,  a  horizontal  resolution  of  60  m  is  possible.  These  resolutions  pre¬ 
suppose  a  sufficiently  intense  laser  pulse  to  overcome  any  noise  problems. 
In  many  current  systems,  it  is  necessary  to  add  several  range  bins  or 
average  several  pulses  together  to  achieve  an  adequate  signal-to-noise 
ratio.  This  decreases  the  resolution  proportionately. 

Lidars  can,  in  principle,  measure  ionospheric  density  at  any 
location  in  view  of  the  satellite.  They  are  not  limited  to  viewing  the 
subsatellite  track,  or  just  the  topside  of  the  ionosphere,  as  are  topside 
sounders.  They  have  the  potential  of  providing  a  three-dimensional  view 
of  the  ionosphere  over  a  large  portion  of  the  earth's  surface.  However, 
the  lidar  system  would  have  to  operate  at  more  than  one  wavelength  to 
compensate  for  the  different  ionic  composition  in  the  various  layers  of 
the  ionosphere.  This  could  require  two  or  more  separate  laser 
transmitters  (but  the  same  receiving  telescope)  if  the  operating 
wavelengths  are  far  different.  Furthermore,  it  may  be  difficult  to  see, 
for  example,  the  lower  side  of  the  F-layer  if  that  layer  is  sufficiently 


dense  to  appear  nearly  opaque  to  the  lidar  pulse.  Higher  lidar  pulse 
power  or  the  use  of  a  different  ionic  transition  and  lidar  wavelength  may 
mitigate  the  problem. 

Lidars  also  have  significant  advantages  over  passive  detection 
systems.  Passive  systems  normally  cannot  provide  electron  density  as  a 
function  of  range  but  can  give  only  the  integrated  or  column  electron 
density  along  the  line  of  sight.  Passive  systems  also  must  rely  on  some 
ambient  excitation  mechanism  which  may  be  variable  or  limited  to  certain 
times  of  the  day.  Lidars  provide  a  constant  excitation  pulse  and  can, 
under  some  circumstances,  operate  night  or  day  with  perhaps  some  degrada¬ 
tion  in  signal-to-noise  ratio  and  resolution  during  the  daytime. 

In  dealing  with  lidars  observing  the  ionosphere,  it  is  first 
necessary  to  determine  the  ionic  species  to  be  observed  and  the 
appropriate  transition  for  that  species.  Figure  8  shows  that  O'*’  is  the 
predominant  species  in  the  F-layer  during  both  day  and  night.  Any  lidar 
instrument  must,  at  least,  be  able  to  measure  this  species  if  it  is  to 
provide  F-layer  electron  density.  Other  species,  notably  NO'*’  and  O2'*’, 
predominate  in  the  lower  regions  of  the  ionosphere  and  must  also  be 
measured.  But  0^  is  both  crucial  and  difficult  to  detect.  Thus,  it  has 
been  considered  to  the  exclusion  of  other  species. 

The  only  transitions  of  0^  appropriate  for  a  lidar  system  are 
the  resonant  or  ground  state  transitions  because  these  transitions  are 
strong  and  rapid.  Strong  transitions  are  necessary  so  that  a  measurable 
portion  of  the  lidar  pulse  is  scattered  back  to  the  detector.  Rapid 
transitions  are  necessary  if  the  lidar  is  to  obtain  range  information. 
Table  7  shows  a  partial  list  of  O'*'  allowed  transitions  with  the  longest 
wavelength  resonant  transitions  circled.  As  previously  mentioned,  these 
all  occur  in  the  far  ultraviolet  region  known  as  the  vacuum  ultraviolet 
(VUV). 


profiles  for  the  midlatitude  ionosphere  at  sunspot 
ng  the  day  and  night  (from  Reference  23). 


allowed  transitions 


To  evaluate  the  potential  of  a  lidar  technique  operating  in  the 
VUV  at  83.3  nm,  it  is  necessary  to  consider  the  scattering  strength  of  the 
O'*”  resonant  transition  relative  to  competing  effects  such  as  total  absorp¬ 
tion  by  0,  O2,  or  N2.  (It  turns  out  that  Rayleigh  and  Mie  scattering  are 
not  competitive  in  the  ionosphere.)  From  standard  spectral  theory^s^  it 
can  be  shown  that  the  atomic  scattering  coefficient  near  an  absorption 
line,  Oy.  is  given  by 


0 


V 


ite^  ^  Y  _ N _ 

me  4ir^  (v-Vq  (y/4ti 


(4-1) 


where  y  is  the  damping  constant,  fgfjs  is  the  absorption  oscillator 
strength,  Vq  is  the  frequency  of  the  absorption  line,  and  N  is  the  den¬ 
sity  of  the  atoms.  The  integrated  scattering  coefficient  near  this  line  is 


^bs^ 


(4-2) 


Equation  4-1  specifies  an  extremely  narrow  line  with  a  width  based  solely 
on  the  lifetime  of  the  upper  state.  In  a  real  situation,  such  as  the 
ionosphere,  thermal  motions  of  the  atoms  create  considerable  Doppler 
broadening  of  spectral  lines.  The  Doppler  line  shape  is  basically  Gaus¬ 
sian  and,  when  combined  with  the  natural  line  shape  given  in  Equation  4-1, 
yields  the  Voigt  function.  This  is  a  difficult  function  to  handle  and 
unnecessary  for  our  purposes.  A  Gaussian  shape  is  adequate  near  the  line 
center  where  a  lidar  system  would  operate.  In  this  case,  the  scattering 
coefficient  can  be  written 

e^  X^Nf 

a  =  /F - 1 —  exp  {-(AX/AXp)^}  (4-3) 

me  2  AAq 
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which  maintains  the  normalization  given  in  Equation  4-2.  The  Doppler 
wavelength,  aXq,  is  given  by 


AXr 


°  Mc2 


(4-4) 


where  M  is  the  mass  of  the  ion,  and  T  is  the  absolute  temperature. 

For  0'*’  at  the  resonant  transition  of  83.3  nm  and  for  a  typical 
ionospheric  temperature  of  1000*K, 


aXq  =  2.83  X  10"^^  cm 
=  .28  pm 

Consequently, 

=  1.23  x  cm^/ion 
N 


(4-5) 


(4-6) 


This  is  a  very  large  cross  section  compared  to  various  total  absorption 
cross  sections  at  the  same  wavelength  for  prominent  neutral  species  in  the 
F-region  (see  Table  8).  O2  has  the  largest  of  the  absorption  cross  sec¬ 
tions  because  83.3  nm  is  close  to  the  first  ionization  edge.  In  contrast, 
the  first  ionization  potential  of  N2  is  only  excited  by  a  shorter  wave¬ 
length  photon. 

Table  9  shows  how  the  resonant  scattering  and  absorption  pro¬ 
cesses  compete  at  various  altitudes  from  300  km  (peak  of  the  F-layer)  to 
200  km  (bottom  side  of  the  F-layer)  to  100  km  (E-layer).  At  300  km  the 
resonant  attenuation  coefficient  is  far  stronger  than  any  of  the  absorp¬ 
tion  coefficients.  These  absorption  coefficients  are  sufficiently  small 
so  that  they  can  be  ignored  here.  At  200  km,  the  resonant  scattering  is 
weaker  and  the  absorption  stronger,  but  the  resonant  scattering  is  still 
predominant.  Below  200  km,  the  O*  density  drops  rapidly  so  that  resonant 
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Table  8.  Cross  sections  at  83.3  nm 


Resonant  Scattering: 


O'*’  —  =  1.23  X  10" 

N 


Total  Absorption: 


0 


^  =  2.6  X  IQ-i 
N 


O2 


—  =  1.5  X  10*1 
N 


N2 


=  2.3  X  10*1 
N 


1^  cm^/ion 


®  cm^/atom 


^  cmVmolecule 


®  cmVmolecule 
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Table  9.  Attenuation  coefficients. 


At  300  km: 

0+  »  loVcm^ 

0  =  7.5  X  10^ 

O2  =  4.4  X  10^ 

N2  =  8.6  X  10^ 

At  200  km: 

0+  =  loVcm^ 

0  =  5.5  X  10^ 

O2  =  2.1  X  108 

N2  =  2.6  X  10^ 

At  100  km: 

O'*'  “  0/cm^ 

0  =  4.1  X  IQi* 

O2  =  2.0  X  10^2 

Nj  =  8.7  X  IOI2 


=  1.23  X  10“^  cm~^ 

oa  =  1.95  X  10-5 

Og  =  6.53  X  10” ' ^ 

Oa  =  1.99  X  10-10 

Oy  =  1.23  X  10"^  cm"i 

Oa  =  1.44  X  10-8 

Oa  =  3.14  X  10-5 

Oa  =  5.98  X  10-5 

Oy  “  0  an~  1 

Oa  =  1.08  X  10-8 

Oa  =  2.99  X  10-8 

Oa  -  2.00  X  10-8 
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scattering  disappears  at  83.3  nm.  In  about  the  same  region,  absorption  by 
the  neutral  species  increases  rapidly  so  that  most  photons  are  absorbed. 
At  100  km  the  photon  intensity  drops  by  1/e  in  .5  km.  Hence,  the  atmos¬ 
phere  appears  black  to  83.3  photons  at  100  km. 

This  absorption  at  100  km  is  useful  because  it  strongly  limits 
background  problems  and  because  it  eliminates  any  radiation  hazard  to 

people.  At  night,  the  background  near  83  nm  should  be  almost  non-existent 
except  perhaps  in  auroral  regions.  During  the  day,  the  background  for  a 
downward  looking  satellite  sensor  should  come  from  resonant  scattering  of 
sunlight  by  O’*".  The  flux  of  sunlight  at  83  nm  is  about  7.4  x  10®  photons/ 
cm^-sec-nm.  A  lidar  of  modest  capabilities  can  compete  with  this  back¬ 
ground  enabling  daytime  as  well  as  nighttime  measurements.  Some  high 

power  laser  systems  that  could  be  put  on  a  satellite  pose  potential  eye 
hazards  if  pointed  at  the  earth.  But  strong  atmospheric  absorption  at  83 
nm  affords  more  than  ample  protection  against  such  hazards  for  the  lidar 
discussed  here. 

When  the  O'*’  density  is  near  or  above  10®/cm®,  a  lidar  looking  at 
the  0^  line  will  have  difficulty  seeing  past  the  peak  of  the  F-layer. 
Table  9  shows  that  I/02  =  8  km  for  [O'*’]  =  10®/cm®.  This  means  that  the 
peak  of  the  F-layer  can  become  optically  thick  to  the  lidar  pulses  and 
block  the  view  of  the  bottom  side.  This  problem  would  be  most  pronounced 
where  the  F-region  is  strongest  (during  solar  maximum  in  the  daytime  near 
the  equator)  and  less  pronounced  to  non-existent  where  the  F-region  is 

weak  (during  solar  minimum  at  night  away  from  the  equator).  Additional 

lidar  pulse  power  can  compensate,  to  some  extent,  for  this  problem. 

Because  no  current  tunable  lidar  system  works  in  VUV  at  83  nm 
and  because  no  lidar  has  been  flown  on  a  spacecraft,  there  is  a  substan¬ 
tial  question  of  the  practicality  of  this  technique.  Current  ground  based 
lidars  are  large  and  heavy,  require  substantial  average  power,  and  are 
tempermental .  However,  lidar  technology  is  rapidly  evolving  so  that  it  is 
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reasonable  to  believe  that  a  practical  lidar  system  similar  to  the  one 
discussed  here  will  be  available  within  the  near  future. 


NASA  has  long  recognized  the  potential  of  space-borne  lidar s  for 
many  types  of  atmospheric  observations  and  now  feels  that  the  technology 
has  advanced  sufficiently  to  propose  an  ambitious  lidar  research  program 
for  the  Space  Shuttle.  Table  10  shows  a  partial  list  of  lidar  experiments 
that  NASA  hopes  to  pursue.  Those  identified  with  arrows  are  particularly 
pertinent  to  the  ionosphere.  It  is  worth  noting  that  the  shortest  wave¬ 
length  under  consideration  is  about  200  nm,  which  is  at  the  boundary  of 
the  VUV,  not  far  from  the  80  nm  necessary  to  work  with  O'*". 

Figures  9  and  10  show  drawings  of  a  lidar  instrument  and 
integration  of  that  instrument  into  the  Space  Shuttle.  Although  the  lidar 
equipment  shown  is  large,  it  is  comparable  to  other  equipment  NASA  intends 
to  fly.  The  lidar  power  requirements,  however,  take  a  disproportionate 
share  of  that  available  to  all  experiments. 

One  of  the  more  interesting  lidar  experiments  planned  for  the 
Space  Shuttle  involves  looking  for  Mg'*'  in  the  E  region  using  an  Mg'*'  reso¬ 
nant  transition.  Table  11  presents  lidar  parameters  which  NASA  expects  to 
achieve  and  which  would  yield  meaningful  data.  Most  of  these  parameters 
are  typical  of  lidars  in  operation  today.  The  factor  limiting  performance 
is  the  pulse  energy  of  .002  joules.  Even  with  this  limitation,  Figure  11 
shows  that  adequate  resolution  is  possible  under  a  variety  of  conditions. 

A  lidar  looking  for  Mg'*'  at  279  nm  has  much  the  same  problem  as  a 
lidar  looking  for  0'*'  at  83  nm.  The  resonant  scattering  cross  section  for 
0'*'  (1.23  X  10"^^  cm^/ion)  is  a  little  smaller  than  for  Mg'*'  (4.27  x  10~^^ 
cm^/ion),  but  the  density  of  O'*"  is  usually  several  orders  of  magnitude 
higher  than  Mg^.  Furthermore,  the  daytime  solar  background  near  the  0'*' 
line  is  more  than  an  order  of  magnitude  smaller  than  near  the  Mg'*'  line. 


50 


10.  Partia]  list  of  lidar  experiments  that  NASA  considers  feasible 
for  the  Space  Shuttle.  Arrows  indicate  those  of  interest  to 
ionospheric  studies  (adapted  from  Reference  22). 


EXPtiUMtNI  CLASS 


WAVELtNCtH  •  NM 


LA.^  HARPWaIU 


SVSli-M4 


CIOLD  TOP  HEIGHTS  A  Sl  RFACt  ALBEDO 
THOPOSPHERK  AEKOKIL  PROFtUS 
SOCT  C LOUD  A  CIRC  IMPOLAR  PART 
LAYER  PROFILE 

STRAT  aEROM>1  BaCRSCAT  PROI  ILl 

aUali  atom  profiles 


lONOSPHERK  MITaL  IO^iS 


CHECMK  AL  RELEASE  OlACNOCTlCs 


STRATOSPHERIC  020WE  PROFILES 


mESOPHERIC  oh  PROFILE 
THE  RMOSFMERIC  SO  PROFILE 


mOPOSPHERT  H:0  PROFILE 


SURFACE  AND  CLOUD  TOP  PRESSURE 
TROPOSPHERIC  PRESSURE  PROFTLi 
TROPOSPHERIC  TEMPERATURE  PROFlU 

TROPOSPHERIC  NO;  PROFILE 


COLUMN  CONTLNT  AND  ROUGH 
PROFILES  OF  following  SPECIES 
FMi  NHidOAvMi  CFM4 
(iC.TRMi.CjH4.  CCI4.I  jHtCI.H.-O  ETC 

STRATOSPHERIC  SPECIES  MEASURE 
MINT  AT  TANCEVT  AlTTrUOt 
C»  NHj  CFM  iiHNOi  CIO  HjHa  CCU 
riONGw  ETC 

(TWO  Satellite  tech  i 


aRRUSKE'WATERDISC 


THERMOSfWERK  M(  Mf  •  acw 
MpO  PROFILE 


SODOM  TEMPERATURE  AND  WINDS 


T6«>«k> 

t’O.’lln 


TMC. 

TAU(.'. 

T«0(:> 


:>«.*>  2»5,:«M50P 


Hd  Vac  FREQUENCY  IMH'BLLD 


NARROW  RAND  n:SABU  DYI 
LAST  iR  pi  MPU)  WITH  Nil  ^a;. 
DOUBLED 

NARROW  band  TtiNABLE  D^  1 
UOtBlt  PUMPED  Wirtl  NJ  T  AC 
DOUBLED 

NARROW  RANOTDNABLt  UVI 
LASER  PI >MPfD  WITH  Nil  TAi. 

DOUBLED 

NARROW  BAND TUKAB'.I  !»( 
LASER  PUMPTP  WITH  Nd  >A1. 
TRIPl  ED 

NARROW  hand  tunable  DYI 
LASER  DOUBLED  PIMPED  WITH 
SU  VaG  DOUBLED 


TWO  NARROW  BANU  TlfNAKt  I 
DYE  LASERS P( MpED  WITH 
Nd  VAC  DOUBLED 


cw4r,»  niscRfriiv  tln'Bli 
WITHlil  TIRCM'VNI  DITECTKin 
TWO  WAS  firNoTHS 


PULSEDCO*  DI.WRETEIY 

riMABti  hTterodym 
OCTEa  TKIN 
TWfiWAVtlF.NCTW> 


■Na  VAO  DOUBLED -POLARIZlL- 


THREI  NARROW  BAKDTUNANlI 
OVr  LASERS  TWO  tXIUBLLD 
AND  PUMPED  WITH  Na  YAG 
TRiriED  AND  ONE  PUMPED  WITH 
Nd  v*r  noURlFD 


NARROW- Band  TI'NaBLI  dm 
lASER  PUMPCO  WITH  Nd  >A( 
OOllBLEIi 

SPtCIAl  VIMRl  E  I  aM  P  I  ON« 
rULSl  >.'»«.''TR>  NARR(>» 
BAND 


SPECIAL  LASER  MODE  L«r  RED 


MNCH  PMI 


SlNI.lt  PM  I 
.MN(;LI  PM’ 


SINGl  t  P 
SINGU  P 
-SISK.LI  P 


UHlLt  PMl  WITH 
POI  AR17MION 
sEPaRaiiov 


•  ■  ;<  MtllR  Kt  I  llv'Fk 
I  OkMIT  A())<  viabII 
S|{  UN|)4M> 

>  .  iiu  RAM.) 


CRYOCINK  INTIGRaTII' 
WITH  LASI  F 
fTWr)DFTirTOR.^< 


ORBIT  ADIl'SI  ABLI 
LASI  R  MOl  M 
NtAR.SIMli  UNEDIS 
laser  PI  LSk*^ 


.SWING  AWAV  HOI  H 
INC.  MIRROR  OR  LASI  R 
transmit  RtCl  IM 
THROU  GH  DtTvnOR  mRI 
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SMROUII 


Figure  9. 


Figure  10. 


Proposed  lidar  instrunent  for  the  Space  Shuttle  (adapted  from 
Reference  22). 


A  lidar  instrument  integrated  into  the  Space  Shuttle  (adapted 
from  Reference  22). 
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NASA  specifications  for  a  Nd:YAG-Dye  laser  system  to 
be  floMn  on  the  Space  Shuttle  to  measure  in  the 
ionosphere-  (Adapted  from  Reference  26). 


Table  11. 

Wavelength 
Laser  Line  Width 
Pulse  Energy 
Pulse  Rate 
Power  Consumption 
Receiver  Area 
Receiver  Bandpass 
Detection  Efficiency 
Receiver  FOV 
Shuttle  Altitude 
Mg'*"  Back  scatter  Strength 
Daytime  Background 


279.6  rrni 
.6  pm 

.002  0  {Limiting  Factor) 

30  Hz 

2  kW 
10“*  cm^ 

1.5  nm 

.02 

.2  mrad 
270  km 

3.4  X  10"^^  cm^  ion"'  sr"' 

3  X  10^  photons  nm" '  sr" '  s" '  cm 
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Vertical  and  horizontal  resolution  tradeoffs  (assuning  a  ±10% 
density  measurement  error)  for  the  proposed  NASA  lidar 
experiment  (from  Reference  Z6). 


Hence,  a  lidar  operating  at  83  nm  with  parameters  similar  to  those  of  the 
proposed  NASA  lidar  will  produce  as  good  or  better  0^  measurements  than 
the  Mg'*"  measurements  anticipated  for  the  NASA  instrument. 


THOMSON  SCATTER  RADARS 


Thomson  or  incoherent  scatter  radars  are  very  useful  for  deter¬ 
mining  electron  density  as  well  as  other  ionospheric  parameters.  However, 
this  type  of  radar  has  only  been  used  for  ground  based  observations  and 
never  on  a  satellite  because  of  the  need  for  high  power  and  large  anten¬ 
nas.  Figure  12  shows  one  ground-based  facility.  With  the  advent  of  the 
Space  Siiuttle  and,  consequently,  the  ability  to  deploy  large,  heavy,  and 
complex  instrumentation  in  space,  such  a  space-borne  radar  becomes  more 
practical.  Power  is  still  a  problem  but  may  be  solved  by  the  deployment 
of  a  small  nuclear  reactor.  Los  Alamos  National  Laboratory  is  reportedly 
designing  such  a  reactor  for  space  application,  and  the  Soviet  Union  is 
known  to  have  deployed  a  reactor  to  power  a  satellite-borne  surveillance 
radar.  (The  nuclear  powered  satellite.  Cosmos  954,  re-entered  the  earth's 
atmosphere  over  northern  Canada  on  the  morning  of  January  24,  1978). 

Conventional  ionospheric  sounding  discussed  in  Section  3  depends 
on  the  reflection  of  radio  waves  by  an  electron  gas.  The  reflection 
process  uses  the  collective  behavior  of  the  electrons  in  a  magnetic  field 
and  is  described  by  the  Appleton-Hartee  equation.  This  is  the  usual  over- 
dense  reflection  which  occurs  when  the  radio  frequency  equals  the  plasma 
frequency,  fp,  given  by  Equation  3-1.  At  frequencies  above  the  plasma 
frequency,  the  radio  waves  will  penetrate  through  the  plasma  essentially 
unimpeded  unless  they  encounter  irregularities  in  plasma  density,  especi¬ 
ally  with  sharp  gradients.  In  this  case,  some  underdense  scattering  will 
occur.  Underdense  echoes  can  sometimes  be  seen  by  topside  sounders  but 
are  usually  weaker  than  overdense  reflections.  For  example,  if  the  plasma 
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Arecibo  Ionospheric  Observatory  v*«1ch  is  often  used  as 
Thomson  scatter  radar  to  probe  the  ionosphere. 


contains  density  irregularities,  aN/N,  that  are  assumed  to  be  spherically 
sytnmetric  and  given  by  a  Gaussian  autocorrelation  function  with  scale 
length  L,  then  it  can  be  shown  that  the  radar  scattering  cross  section  per 
unit  volume  is  given  by^^. 28 


s  -  ^  l-l'  f'  (2.)^^  L’  exp[.8.(i)'l  .  (5-1) 

N  ^  X 

This  cross  section  depends  strongly  on  the  value  of  L/x  appearing 
in  the  exponential  term  and  will  become  very  small  for  x  «  L,  Typical 
ionospheric  irregularities  associated  with  spread-F  may  have  aN/N  «  .01 
with  L  of  the  order  of  10  m.  In  this  case,  underdense  reflections  will 
virtually  disappear  for  x  <  1  m  unless  still  smaller  density  fluctuations 
exist.  Because  of  the  random  thermal  motions  of  the  electrons  and  ions, 
these  smaller  fluctuations  do  exist.  The  scattering  caused  by  these  fluc¬ 
tuations  is  usually  referred  to  as  "Thomson”  or  "incoherent"  scatter. 
Both  of  these  terms  are  in  common  use  but  both  are  technically  incorrect 
because  of  the  collective  effects  that  occur. 


After  J.  J.  Thomson  discovered  the  electron,  he  showed  that  a 
radio  wave  falling  on  a  free  electron  will  cause  it  to  oscillate  and 
thereby  to  radiate  a  scattered  wave  in  all  directions. 29  The  power 
reradiated  in  this  way  is  equal  to  the  power  falling  on  an  area  given  by 


8tt 


me  2 


(5-2) 


where  rg  =  e^/mc2  =  2.82  x  10" cm  and  is  known  as  the  classical  radius 
of  the  electron.  oj  is  the  total  isotropic  scattering  cross  section 
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of  a  single  electron  and  is  called  the  Thomson  cross  section.  The 
differential  scattering  cross  section  can  be  written  as 

°e  ^ 

where  iii  is  the  angle  between  the  direction  of  the  incident  electric  field 
and  the  observing  direction  ( tp  =  ir/2  for  backscatter) .  Note  that  both 
Equations  5-2  and  5-3  are  independent  of  frequency. 

Because  the  electrons  in  any  volume,  V,  have  random  thermal 
motions,  they  will  scatter  radio  waves  in  a  random  fashion  imparting  a 
random  phase  to  them.  This  gives  rise  to  the  term  "incoherent"  scatter. 
A  receiver  looking  for  these  scattered  waves  will  receive  a  signal  power, 
Pf^,  representing  the  sum  of  signals  with  random  phases  and,  therefore, 
proportional  to  the  density  of  electrons  in  the  volume  illuminated.  If 

this  volume  is  at  a  range,  r,  from  a  transmitter  with  pulse  power,  Pj, 

the  power  returned  per  unit  area  is  given  by  the  familiar  radar  equation 

’’r  "  ^  o^NV/[4n  r2]2  .  (5-4) 

Hence,  by  measuring  the  returned  power,  the  electron  density  N  may  be 
determined.  The  volume,  V,  must  be  determined  from  the  geometry  of  the 

experiment.  The  Thomson  scattering  cross  section  of  an  electron,  a,  is 

very  small  (>=>10"2‘^  cm^)  so  that  the  power  returned  from  the  ionosphere  is 
also  very  small. 

The  electrons  that  scatter  incident  radio  waves  move  with  veloc¬ 
ities  which  depend  on  their  temperature  causing  a  Doppler  shift  in  the 
scattered  wave.  The  Doppler  shift  appears  as  a  substantial  frequency 
broadening  in  the  returned  pulse.  The  scattered  spectrum  would  be 


60 


expected  to  have  a  Gaussian  shape  with  a  full  width  at  half  maximum  of 
1.42  Afg  where  Afg  is  defined  as  the  Doppler  shift  of  an  electron 
approaching  the  receiver  at  a  mean  thermal  speed.  This  is  given  by 

Afg  =  1  (8kTg/mji/2  (Hz)  (5-5) 

where  k  is  the  Boltzmann  constant,  and  Tg  is  the  absolute  electron 
temperature. 

W.  E.  Gordon  was  the  first  to  demonstrate  that  this  Thomson 
scattering  might  be  observed  by  radar  techniques  and  suggested  the  use  of 
a  wavelength  of  x  =  1.5  m.^o  When  Tg  ■»  1600*K,  1.42  Afg  «  500  kHz. 
Fortunately,  this  very  large  bandwidth  did  not  materialize  when  actual 
experiments  were  conducted.  K.  L.  Bowles  observed  a  scattering  similar  to 
Thomson  scattering,  but  confined  to  a  much  narrower  bandwidth.  This  was 
correctly  attributed  to  an  electrostatic  coupling  between  the  electrons 
and  ions. 

When  the  incident  wavelength  is  much  larger  than  the  Debye 
screening  distance  given  by22 

D  =  (eokTg/4„Ne2)i/2  (5-6) 

where  gq  is  the  permittivity  of  free  space,  and  e  is  the  charge  on  the 
electron,  the  electrons  show  a  collective  behavior,  and  the  scattering 
becomes  partially  coherent.  Because  of  the  coupling  of  the  electrons  to 
the  ions  and  the  much  slower  thermal  velocity  of  the  ions,  the  scattering 
bandwidth  is  much  smaller  than  originally  predicted  by  Reference  30.  The 
Debye  distance  is  typically  1  cm  or  less  in  the  ionosphere  compared  with  a 
typical  radar  wavelength  near  1  m.  Hence,  partially  coherent  scattering 
should  occur. 


Pure  incoherent  Thomson  scattering  is  still  possible  when  x  <  D, 
but  it  is  exceedingly  impractical  because  of  the  broad  scattering 
bandwidth  that  results.  Such  scattering  has  been  observed  in  laboratory 
experiments  using  a  laser  instead  of  a  radar.  Use  of  a  laser  for  Thomson 
scattering  is  possible  with  the  ionosphere  but,  probably,  impractical. 
Unlike  the  1 idar  system  proposed  in  Section  4,  this  technique  would  have 
the  virtue  of  yielding  a  direct  measurement  of  electron  density,  not  just 
ion  density. 

In  the  usual  case  where  x»  D,  the  collective  behavior  of  the 
electrons  and  the  coupling  of  electrons  to  ions  offers  the  possibility  of 
obtaining  additional  information  from  the  scattered  signal.  Not  only  can 
electron  density  be  determined,  but  also  ion  composition  if  the  ratio  of 
electron  to  ion  temperature  is  known.  Conversely,  if  only  a  single  ion 
species  is  present,  the  ratio  of  electron  to  ion  temperature  can  be 
determined. 32^ 33  Also,  it  is  possible  to  use  the  Faraday  rotation  of  the 
scattered  signal  to  determine  the  column  electron  density.  (This  latter 
technique  is  discussed  in  a  different  context  in  the  next  section.) 

When  X  >  0,  the  scattering  cross  section  given  by  Equation  5-3 
must  be  modified  to  give,^*^ 

4„(r  sin^)2 

o  =  - -  (5-7) 

(Ua^)(UT^/T.+«2) 


where 


o  -  4it  D/x  <  1  .  (5-8) 

For  typical  ionospheric  parameters,  the  denominator  of  Equation  5-7  is 
about  3.  This  means  that  the  scattering  cross  section  in  this  case  is 
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smaller  than  the  Thomson  cross  section  by  about  a  factor  of  3.  Because  of 
the  much  narrower  bandwidth  of  the  scattered  signal,  however,  the  scatter¬ 
ing  represented  by  Equation  5-7  is  far  easier  to  detect  than  pure  Thomson 
scattering.  But  the  detection  of  any  return  is  still  difficult  because  of 
the  very  small  cross  section  involved  (typically  3  x  cm^). 

Figure  13,  taken  from  a  study  of  Thomson  scatter  radars  for  the 
Space  Shuttle, 3**  shows  anticipated  noise  power  as  a  function  of  frequency, 
pulse  width,  and  whether  the  antenna  is  pointing  toward  or  away  from  the 
earth.  This  figure  shows  that  frequencies  of  several  hundred  MHz  exhibit 
the  lowest  noise  power  for  an  antenna  pointed  away  from  the  earth.  For  an 
antenna  pointed  toward  the  earth,  the  best  choice  would  be  frequencies 
below  100  MHz.  Figure  14  is  similar  to  Figure  13  and  shows  anticipated 
signal-to-noise  ratio  versus  frequency.  The  conclusion  as  to  optimum 
frequency  remains  the  sane. 

Another  consideration  is  the  required  power-aperture  produce. 
To  simplify  consideration  of  this,  the  radar  wavelength  is  assumed  to  be  \ 
»  .5  m  (600  MHz)  which  is  near  the  optimum  indicated  by  Figures  13  and 
14.  Furthermore,  it  is  assumed  that  a  signal-to-noise  ratio  of  0.1  is 
required  at  a  range  of  100  km  with  an  electron  density  of  10^  electrons/ 
cm^.  In  this  case,  it  can  be  shown  that  the  radar  must  have  a  power- 
aperture  produce  of  at  least  10^  watts-m^,  and  that  a  larger  one  would  be 
preferable. 3'*  Figure  15  illustrates  the  tradeoffs  between  power-aperture, 
pulse  width,  and  antenna  pointing  direction. 

To  achieve  such  a  large  power -aperture  produce  requires  both  a 
high  pulse  power  and  a  large  antenna.  The  largest  antenna  operating  near 
600  MHz  proposed  for  the  Space  Shuttle  is  183  m  in  diameter.  This  has  an 
effective  area  of  about  1.3  x  10**  m^.  Assuming  a  power  aperture  product 
of  10®  watts-m^,  implies  tha  the  peak  pulse  power  must  be  10  kW.  This  is 
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Figure  13.  Anticipated  noise  poi«er  versus  frequency  for  a  Thomson  scatter 
radar  In  orbit  and  pointing  either  toward  or  away  from  the 
earth,  t  is  the  transmitted  pulse  width.  Shorter  pulse 
widths  require  a  larger  receiver  bandwidth  which  in  turn 
admits  more  noise  at  a  given  frequency  (from  Reference  34). 


Power-aperture  product  required  in  a  typical  case  (from 
Reference  34). 


well  within  reason  if  the  duty  cycle  is  not  too  large  because  the  average 
power  requirement  would  be  considerably  less.  More  power  may,  in  fact,  be 
possible.  At  a  peak  power  of  one  megawatt,  however,  antenna  breakdown 
could  be  a  problem. 

To  illustrate  the  size  of  a  183  meter  antenna.  Figure  16  shows 
one  unfurled  from  the  Space  Shuttle.  One  difficulty  with  such  a  large 
antenna  is  the  drag  that  it  encounters  in  orbit.  Figure  17  shows  the 
expected  altitude  decrease  with  time  for  various  size  antennas  starting  at 
a  300  km  altitude.  It  is  clear  that  a  600  ft  (>=183  m)  antenna  would  not 
remain  in  a  300  km  orbit  very  long,  but  would  remain  much  longer  in  a 
higher  orbit. 

Placing  the  radar  in  a  higher  orbit,  however,  increases  the 
range  to  the  ionosphere  and  decreases  the  signal-to-noise  ratio.  Figure 
18  shows  this  decrease  with  range.  Note  that  the  configurations  are  for  a 
10  kW  pulse  power.  A  higher  pulse  power  could  compensate  for  this 
decrease  of  signal-to-noise  ratio  with  distance. 

Although  it  might  seem  that  higher  pulse  power  is  always  desir¬ 
able,  this  may  not  be  the  case  for  a  spacecraft  where  the  average  power  is 
limited.  In  this  case,  it  is  necessary  to  consider  the  trade-offs  between 
pulse  power  and  duty  cycle.  Higher  pulse  power  increases  the  signal-to- 
noise  ratio  while  a  higher  duty  cycle  allows  more  measurements  and  a 
better  determination  of  average  ionospheric  characteristics.  Regardless 
of  the  signal-to-noise  ratio,  it  is  necessary  to  do  some  averaging  of  the 
Thomson  scatter  signal  to  accurately  measure  the  mean  characteristics  of 
the  ionosphere. 

It  can  be  shown  that  the  standard  deviation,  o,  of  the  Thomson 
scatter  measurement  varies  approximately  as 
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hypothetical  radar 

(from  Reference  34).  unfurled  from  the  Space 


where  S/N  is  the  signai-to-noise  ratio,  C  is  a  constant  near  1  and  n  is 
the  number  of  samples  used  in  the  average.  When  S/N  is  small,  a  varies  in 
proportion  to  S/N.  As  S/N  approaches  1,  further  increases  have  little 
effect  on  a.  In  this  case,  it  is  better  to  increase  the  number  of 
samples,  n,  in  order  to  improve  o.  In  practice,  it  is  necessary  to  add 
many  pulses  together  to  minimize  a.  Adding  many  pulses  together  leads  to 
a  relatively  poor  horizontal  resolution  because  a  downward  looking  Thomson 
scatter  radar  at  a  Space  Shuttle  altitude  of  about  300  km  sweeps  through  a 
large  region  of  the  ionosphere  in  a  short  time.  Depending  on  the  design 
of  the  radar,  horizontal  resolution  could  vary  from  50  to  more  than  800 
km. 

Vertical  resolution,  in  contrast,  is  a  function  of  pulse  length 
and  would  vary  from  3  km  to  30  km  for  typical  radars.  This  is  similar  to 
that  of  topside  sounders  while  the  horizontal  resolution  is  much  worse 
than  topside  sounders. 

Another  serious  problem  is  clutter  caused  by  reflections  from 
the  earth's  surface,  the  auroral  region,  or  anything  other  than 
ionospheric  electrons.  Figure  19  shows  an  example  of  a  signal  return  for 
a  single  transmitted  pulse.  The  ionospheric  returns  are  small  compared 
with  the  ground  clutter.  The  ground  clutter  can,  however,  be  range-gated 
out,  but  doing  so  limits  the  pulse  repetition  frequency.  F-region  auroral 
clutter  may  be  similarly  strong  and  not  possible  to  eliminate  by  range- 
gating.  It  is  a  substantial  unknown.  Because  clutter  is  often  strong,  it 
may  be  a  problem  even  for  a  directional  antenna  by  feeding  into  side 
lobes. 


Figure  19.  Anticipated  signal  return  for  a  radar  looking  toward  the  earth 
(from  Reference  34). 


In  sunmary,  Thomson  scatter  radars  in  space  appear  extremely 
cumbersome  although  perhaps  feasible.  They  require  a  large  antenna  and 
high  pulse  power.  Vertical  resolutions  of  3  to  30  km  are  possible,  while 
horizontal  resolutions  of  50  to  800  km  may  be  possible.  Horizontal 
resolution  depends  on  the  duty  cycle  and  average  power  available  from  the 
radar.  Ground  clutter  limits  the  pulse  repetition  frequency.  Auroral 
clutter  is  potentially  a  serious  problem  which  would  render  the  radar  use¬ 
less  in  the  auroral  region. 


SECTION  6 

FARADAY  AND  PHASE  METHODS 


The  phase  of  radio  frequency  signals  is  significantly  modified 
by  propagtion  through  the  ionosphere.  These  phase  changes  as  measured  by 
phase  itself,  rate  of  change  of  phase  (Doppler  shift)  or  difference  in  X- 
and  0-  mode  phase  (Faraday  rotation)  are  directly  related  to  properties  of 
the  electron  density  integrated  along  the  signal  path.  In  this  chapter, 
the  potential  for  using  phase  measurements  in  conjuction  with  other  infor¬ 
mation  to  find  local  electron  densities  along  a  signal  path  is  reviewed. 
In  Section  6.1,  techniques  for  radio  frequencies  well  above  the  peak  iono¬ 
spheric  plasma  frequency  are  considered.  In  Section  6.2,  techniques 
applicable  when  the  radio  frequencies  are  near  the  plasma  frequency  are 
explored.  Finally,  in  Section  6.3,  the  utility  of  these  methods  are 
assessed.  It  is  established  that  these  passive  phase  methods,  which 
generally  operate  between  a  satellite  and  a  ground  station,  are  signifi¬ 
cantly  limited  in  usefulness.  Only  the  region  along  the  1 ine-of-sight 
between  transmitter  and  receiver  is  sampled,  and  the  phase  data  obtained 
is  directly  related  to  the  integrated  electron  density  along  the  path. 
This  data  can  be  used  in  conjunction  with  that  gained  by  other  techniques 
to  deduce  features  of  the  electron  density  profile  but  cannot  stand  alone 
for  this  purpose.  Still,  phase  measurements  are  relatively  easy  to  take 
and  interpret  so  that  they  can  profitably  be  used  in  conjunction  with 
other  monitoring  techniques  as  a  routine  source  of  information  acting  as  a 
constraint  on  other  types  of  measurements. 
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6.1  FIRST  ORDER  METHODS 

6.1.1  Doppler  Shift 

If  a  radio  signal  in  the  ionosphere  has  a  transmitted  frequency, 
f,  much  greater  than  the  plasma  frequency 

fp  =  (Ngr^cViT)i/2  (6_i) 

(where  Ng  is  electron  density,  rg  is  the  classical  electron  radius 
(2.82x10"^^  cm),  and  c  is  the  speed  of  light),  then  the  index  of  refrac¬ 
tion  of  the  ionosphere  is  to  first  order 

M  =  1  -  1/2  (fp/f)2  .  (6-2) 

The  phase  shift  imparted  to  a  signal  passing  through  the  iono¬ 
sphere  is  thus 


Acf)  =  /  (1-u)  ds  (radians) 

"  R 


(6-3) 


where  the  integration  is  along  the  ray  path  from  receiver  to  transmitter. 
Therefore,  to  first  order  (for  which  the  path  is  a  straight  line) 


r  c  T 

A-t.  =  ^  /  Ngds 

^  R 


(6-4) 


If  X  is  the  zenith  angle  of  the  wave  vector  this  can  be  written  as 


r  c  Zj 

= -  /  N  sec  X  dz 

f  ^  e 
2r 


(6-5) 
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where  and  z-j-  are  the  heights  of  the  receiver  and  transmitter. 

Except  for  the  sec  x  factor,  the  integral  in  Equation  6-5  yields  the  total 
electron  content  (TEC)  of  the  ionosphere  along  th'  path.  However,  unless 
the  path  is  vertical,  x  is  not  constant,  so  that  the  TEC  can  be  found 
precisely  from  a  measurement  of  only  if  the  shape  of  the  electron 

density  profile  is  known.  Direct  measurements  of  A((i  are  difficult  so 
often  its  rate  of  change  is  measured.  dA4>/dt  is  the  Doppler  shift  which 
is  nonzero  only  if  there  is  some  relative  motion  between  the  satellite  and 
ground  station  or  if  the  ionosphere  is  not  stationary.  Given  Doppler 
shift  data  in  conjunction  with  known  satellite  motion,  it  is  again 
possible  to  evaluate  the  slant  TEC  as  in  Equation  6-5.  However,  any 
horizontal  gradients  in  Ng  will  influence  the  Doppler  shift  so  these 
gradients  introduce  uncertainties  in  the  slant  TEC.  If  the  ionosphere  is 
stationary  and  has  no  horizontal  gradients  the  Doppler  shift  is 


f  =  (seTT)  .  /T  N  dz  +  i  R  (6-6) 

0  f  dt  e  c 

where  R  is  the  component  of  the  satellite  velocity  along  the  path  and 
sec  X  ■'s  the  mean  value  of  sec  x  weighted  by  Ng.  (Note  that,  at  the 
point  of  closest  approach,  fp  will  vanish  except  for  the  effects  of 
horizontal  gradients  or  nonstationarity. )  Again  the  TEC  cannot  be 
obtained  without  knowledge  of  the  shape  of  the  density  profile. 

6.1.2  Faraday  Rotation 


In  the  absence  of  collisions,  the  index  of  refraction  of  an 
electromagnetic  wave  in  a  magnetized  plasma  is^® 


y  -  [  1  - 


2X(1-X) 


2(1-X)  -  Y^sin^e  ±  [Y‘’sin‘*e+4(1-X) ^Y^cos^e] 


1/2 


1/2 

1  (6-7) 
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where 


X  “  (f/f)" 

Y  •  f  /f 

e 

f  =  (electron  gyrofrequency) 

e  2Trin^, 

and  G  is  the  angle  of  the  wave  vector  with  respect  to  the  magnetic  field 
The  +  is  for  the  "ordinary”  wave  (0)  and  -  is  for  the  "extraordinary" 
wave  (X),  If  the  frequency  is  well  above  ^the  plasma  frequency,  and 
the  propagation  is  not  close  to  transverse  to  B,  then  the  quasi-longitudi¬ 
nal  approximation  applies.  This  gives  to  first  order  in  X 

«  1  -  1/2  X  ±  1/2  X  Y  cose  (6-8) 

To  this  same  order  of  approximation,  the  X  and  0  waves  are  circularly 
polarized  in  opposite  senses^^  and,  as  a  result,  a  plane  polarized  wave 
will  have  its  plane  of  polarization  rotate  as  it  propagates  if  the  indices 
of  refraction  vi+  and  p.  differ.  From  Equation  6-8  it  is  seen  that  this 
Faraday  rotation  is  given  by 

f  T 

«  =  —  /  (m.  -  u  )  ds  (radians)  (6-9) 

R 

which  to  lowest  order  is 

er  T 

9.  =  — /  N  B  cose  ds  (6-10) 

2nmf^  p 


If  X  is  the  zenith  angle  of  the  wave  vector  this  becomes 


(6-11) 


er  T 

n  =  /  *^6  ^  dz 

If"!  Zr 

Hence,  measurement  of  ci  would  yield  the  TEC  if  B,  e  and  x  constant 

along  the  path.  Since  this  is  generally  not  true,  a  Faraday  rotation 
measurement  will  yield  the  TEC  only  if  the  shape  of  the  electron  density 
is  known  as  is  the  earth's  magnetic  field.  As  with  the  phase,  it  is  often 
convenient  to  measure  the  rate  of  change  of  the  Faraday  rotation  which 
gives 


6r  2 

^  =  — 2—  1-  [  N  B  cose  secx  dz  ] 

dt  2nmf^  dt  Zr 


(6-12) 


which  contains  contributions  from  satellite  motion  and  horizontal  gradi¬ 
ents  in  Ng  as  well  as  o  and  x- 

6.1.3  Application  to  Properties  of  Electron  Density  Profile 

6. 1.3.1  TEC  Evaluation.  Many  workers  have  used  both  Faraday  and  phase 

measurements  to  produce  estimates  of  the  TEC  using  sufficiently  high  fre¬ 
quencies  and  proper  geometry  to  allow  the  first  order  quasi-longitudinal 
approximation  to  apply.  Present  techniques ^7  allow  accuracy  of  a  few  per¬ 
cent  in  the  raw  TEC  measurements  at  <.150  MHz.  However,  since  n  and 
measurements  depend  upon  the  shape  of  the  profile  (as  well  as  magnetic 
field  for  n),  this  shape  is  required  in  order  to  obtain  a  unique  TEC. 
Application  of  a  variety  of  techniques  to  this  question  have  been  consid¬ 
ered  by  Millman^®.  Using  a  realistic  model  for  the  ionosphere,  including 
horizontal  gradients,  he  has  found  the  sensitivity  of  the  TEC  estimates  to 
the  choice  of  mean  height  of  the  ionosphere.  He  finds  that  a  variation  of 
■50  km  in  the  mean  ionospheric  height  produces  changes  of  ■lO  percent  in 
TEC  estimates.  This  sensitivity  is  a  function  of  the  geometry  as  well  as 
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the  technique  used.  Still,  finding  a  TEC  from  measurements  of  n  or  A4) 
requires  use  of  an  ionospheric  profile,  and  uncertainties  in  this  profile 
will  result  in  uncertainties  of  40  percent  in  the  resulting  TEC.  Gener¬ 
ally,  a  given  measurement  serves  to  scale  a  presumed  shape  of  density 
profile. 

6. 1.3. 2  Atmospheric  Scale  Height  Evaluation.  Potts^^  has  noted  that 
given  a  measurement  of  dn/dt  and  a  value  of  the  peak  electron  density 
obtained  by  a  separate  ionosonde  measurement,  it  is  possible  to  produce 
some  characteristic  of  the  shape  of  the  electron  density  profile.  Since 
the  peak  in  density  contributes  most  strongly  to  the  Faraday  rotation, 
Potts  has  attempted  to  evaluate  what  is  effectively  the  thickness  of  the 
F-region  peak  by  using  an  a-Chapman  profile  whose  parameters  are  peak 
height,  h^,  and  neutral  scale  height,  H.  Potts  found  that  the  value  of 
Q  is  rather  insensitive  to  variations  in  h^  so  that  H  can  be  found  from 
the  Faraday  rotation  measurement  when  used  in  conjunction  with  a  knowledge 
of  the  peak  electron  density.  Still,  this  work  showed  substantial  sources 
of  error  in  horizontal  gradients,  ionospheric  irregularities  and  devia¬ 
tions  from  the  ideal  Chapman  profile.  When  the  peak  electron  density  was 
fixed  by  an  ionosonde  600  km  distant,  it  was  found  that  H  had  an  uncer¬ 
tainty  of  roughly  ±  25  percent. 

6. 1.3. 3  Profile  Parameters  by  Multiple  Measurements.  Davies,  Fritz  and 
Gray 3''  have  reported  simultaneous  measurements  of  both  modulation  phase, 
({,i^  (similar  to  group  path)  and  Faraday  rotation  from  140  MHz  transmis¬ 
sions  of  ATS-6.  These  measurements  are  significantly  different  for  a 
geostationary  satellite  because  the  earth’s  magnetic  field  is  greatly 
reduced  at  the  orbital  radius  so  that  Faraday  rotation  is  largbly  unaffec¬ 
ted  by  plasmaspheric  electron  density  (i.e.,  above  2000  km).  In  order  to 
translate  q  and  ^  measurements  into  electron  content,  use  was  made  of 
monthly  median  values  of  peak  Ng  as  obtained  by  nearby  ionosondes  in 
conjunction  with  an  ionospheric  model  devised  earlier  by  Titheridge*^®. 


Attempts  to  use  concurrent  ionosonde  data  failed  due  to  substantial  and 
variable  horizontal  gradients  over  the  several  hundred  kilometers  between 
the  ionosondes  and  penetration  point  for  the  rays.  No  attempt  was  made  to 
obtain  TEC  values  from  individual  measurements.  Exclusively  monthly 
median  parameters  were  established.  The  model  electron  density  profile  was 
parametrized  with  five  quantities:  height  of  maximum  Ng,  ratio  of  [H'*"  ] 
to  [O'*’]  at  the  peak,  ion  temperature  at  the  peak,  electron  temperature 
at  the  peak  and  the  electron  temperature  gradient  above  the  peak.  These 
parameters  (only  one  of  which  directly  relates  to  the  profile)  were  varied 
within  limits  imposed  by  physical  considerations  in  order  to  obtain  a  fit 
to  the  measured  values  of  and  4^  using  a  trial  and  error  procedure. 
It  appears  that  with  this  many  parameters  using  a  mathematical  best  fit 
would  not  be  suitable.  Using  this  analysis,  Davies^^  finds  that 

uncertainties  of  +3  percent  in  n  and  4^  can  produce  monthly  median 
ionospheric  profiles  which  are  generally  consistent  with  Titheridge's 
model  for  physically  reasonable  parameters.  Normally,  only  a  narrow  range 
of  any  one  parameter  will  allow  choices  of  the  other  parameters  which  are 
reasonable  while  still  fitting  the  observed  median  si  and  4^^.  However, 
the  authors  do  not  give  enough  quantatitive  details  to  measure  the 
possible  variation  in  the  profile  itself. 

When  these  same  techniques  were  applied  to  ATS-6  measurements  in 
India  by  Davies,  Donnally,  Grubb,  Rao,  Rastogi,  Deshpande,  Chandra,  Vats 
and  Sethia**^,  much  less  satisfactory  results  were  obtained  since  the 
propagation  path  crossed  the  magnetic  equator.  As  a  result,  the  Faraday 
rotations  were  greatly  reduced  with  a  corresponding  loss  in  accuracy  to 
±15  percent.  Consequently,  fitting  the  data  with  choices  of  ionospheric 
profiles  was  not  possible  in  this  case  even  though  additional  data  from  40 
MHz  was  used. 
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6.2 


HIGHER  ORDER  CORRECTIONS 


In  Section  A,  only  the  high-frequency  approximation  was  consid¬ 
ered,  that  is  f  »  fp.  Under  this  assumption  the  index  of  refraction  of 
Equation  7  can  be  expanded  to  lowest  order  with  adequate  accuracy.  How¬ 
ever,  if  frequencies,  f,  comparable  with  fp  are  used  for  transiono- 
spheric  propagation,  the  first  order  expansion  is  not  generally  adequate. 
The  appropriate  expression  for  p  under  quasi-longitudinal  (QL)  conditions 
is 


V 


[  1 


X _ t1/2 

1  ±  Ycose 


(6-13) 


which  requires  (QL  condition) 

1  -  X  »  1/2  Y  sino  tane 


(6-14) 


As  f*fp  so  that  X->1  the  QL  applies  only  near  e  =  0.  The  expression  for 
u  can  be  expanded  giving,  for  Y  <  1, 


=  1  -  1/2  X  -  1/8  X2  +  ... 
±Y(l/2  X  +  1/4  X2  +  ...  ) 


(6-15) 


+  higher  order  terms  in  Y. 


Thus,  for  Faraday  rotation 


9. 


p_  )ds 


=  /  (XY  +  1/2  X^Y  +  ...)  ds 

R 


(6-16) 
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L  /  ds 

o  * 
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f  f2  +  1_ 
e  p  ^4 


R 


Apart  from  constant  factors,  this  can  be  written  to  second  order  as 


n  =  —  /  N  B  cose  ds 
f2  O  " 


T 

/  B  cose  ds 
•'  e 

R 


(6-17) 


This  suggests  that  a  measurement  of  n  may  allow  an  estimate  of  /N^s  to  be 

made  if  a  high  frequency  reference  value  of  k  is  obtained  to  fix  /Ngds. 

This  results  from  the  difference  in  frequency  dependence  for  jN^s.  Yeh'*^ 

has  detailed  the  application  of  this  technique  to  establishing  JN^s,  but 

e 


it  was  not  applied  to  any  data. 


6.2.1  Application  to  Data 

When  the  high  frequency  approximation  is  not  adequate,  it  is 
necessary  to  consider  both  corrections  to  the  index  of  refraction  and  to 
any  deviations  from  a  straight  line  ray  path.  In  early  work  with  20  MHz 
and  40  MHz  signals,  this  was  regarded  as  a  substantial  inconvenience  and 
source  of  uncertainty.  Garriott**^  devised  a  conversion  factor  for  these 
refractive  effects  in  order  to  more  accurately  find  the  TEC.  He  found  that 
if  the  transmission  frequency  was  only  twice  the  peak  plasma  frequency  a 
ten  percent  correction  was  required  in  converting  n  to  TEC  using  the  first 
order  formulation.  Garriott  used  analytic  techniques  requiring  simple 
models  of  the  ionospheric  profile  to  estimate  second  order  corrections. 
This  analytic  approach  has  been  extended  by  Ross*^*^  who  established  the 
expansions  to  higher  orders  suggesting  a  general  formalism  for  incorporat¬ 
ing  corrections.  Potts  used  a  different  philosophy  in  accounting  for 
refractive  effects  at  20  MHz.  Instead  of  expanding  u,  he  used  the  full 
expressions  in  a  ray  tracing  procedure  which  required  no  approximations. 
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The  full  effects  of  X  and  0  path  splitting  were  not  included,  but  Faraday 
rotation  was  evaluated  along  a  single  curved  ray  path.  Such  a  procedure 
excludes  the  possibility  of  directly  connecting  and  TEC.  Potts  found  a 
profile  by  fitting  it  to  u  so  that  the  TEC  could  have  been  deduced 
although  it  was  not  done. 

Davies,  et  al.**^  have  reported  on  corrections  applicable  to  the 
40  MHz  phase  and  Faraday  measurements  from  the  ATS-6  satellite.  Using 
typical  ionospheric  shapes  in  conjunction  with  ray  tracing,  they  express 
an  enpirical  relation  between  ()>„,  and  slant  TEC  as 

Nt  ""  (6-18) 

where  u  typically  produces  a  3  percent  correction  at  40  MHz.  Similar 
results  were  obtained  for  i2.  This  empirical  relation  is  applicable  only 
at  this  frequency  and  for  the  ATS-6  to  Boulder  geometry.  No  effort  was 
made  to  estimate  /N^ds  using  this  technique,  since  the  estimated  uncer¬ 
tainty  in  <j)ni  is  also  about  3  percent. 

6.3  POTENTIAL  VALUE  OF  FARADAY  AND  PHASE  MEASUREMENTS 

Present  methods  allow  measurements  of  the  effects  of  the  iono¬ 
sphere  on  HF  and  VHF  signals  to  be  used  to  find  TEC  to  an  accuracy  of  a 
few  percent.  The  combination  of  Faraday  rotation  and  phase  measurements 
can  be  used  with  a  knowledge  of  peak  Ng  to  restrict  the  electron  density 
profile  which  could  produce  those  measurements.  Still,  this  method  has 
not  been  applied  to  find  the  profile  at  a  specific  time. 

fi.j.i  Use  of  Frequencies  Near  the  Plasma  Frequency. 

■.fun  ' rjuMonospher  1C  signals  are  propagated  at  frequencies,  f, 
..iins  are  not  straight,  and  the  X  and  0  paths  are 


different.  Both  the  Faraday  rotation  and  phase  effects  can  be  expressed 
as  a  series  of  terms  of  the  general  form 


A  /  N^ds 
J  e 


+  ^  /  ds  +  . . . 
f2  ^  e 


(6-19) 


so  that  at  lower  frequencies  there  are  significant  contributions  from  the 
integral  of  higher  powers  of  .  These  higher  powers  give  their  major 
contributions  around  the  peak  so  that  use  of  lower  frequencies  chiefly 
provides  additional  information  about  the  peak  of  the  electron  density 
profile. 


Two  approaches  to  such  a  low  frequency  method  might  be  con¬ 
sidered.  Either  an  attempt  can  be  made  to  better  define  the  profile  near 
the  peak  or  to  find  /N^ds  from  the  measurements.  The  latter  idea  may 
prove  difficult,  because  higher  order  corrections  come  not  only  from  inte¬ 
grals  of  higher  powers  of  Ng  but  also  from  the  increasingly  distorted 
ray  paths  when  frequencies  approach  the  plasma  frequency.  As  a  result,  it 
Is  difficult  to  translate  changes  in  p  or  ()>  into  the  corresponding  inte¬ 
grals  over  the  straight  path  implied  by  /N^s.  Furthermore,  it  will 
always  be  necessary  to  truncate  the  series  leading  to  substantial  uncer¬ 
tainties.  For  example,  if  the  overall  uncertainty  of  the  measured  u  is  ±5 
percent  then  gives  a  10  percent  second  order  correction 

resulting  in  a  ±50  percent  uncertainty  in  the  estimate  for  /N^s. 

It  appears  to  be  better  to  avoid  trying  to  measure  higher  order 
analogues  of  the  TEC  by  directly  using  low  frequency  measurements  to 
refine  the  region  of  the  profile  near  the  peak  in  electron  density  since 
this  is  in  fact  the  quantity  of  ultimate  interest.  This  will  avoid  the 
rather  artificial  problems  associated  with  the  expansion  and  truncation  of 
the  expressions  for  the  path  integrals  in  the  presence  of  complex  paths 
which  themselves  depend  upon  the  ionospheric  profile.  Techniques  for  the 


direct  fitting  of  electron  density  profiles  from  relatively  low  frequency 
measurements  of  Faraday  rotation  and  phase  will  now  be  considered. 


6.3.2  Application  of  Lower  Frequency  Methods 

It  has  already  been  demonstrated  by  Davies that  measurements 
of  >2  and  (j),,!  to  an  accuracy  of  a  few  percent  at  140  MHz  can  be  used,  in 
conjunction  with  peak  density  and  a  physical  model  of  the  ionosphere,  to 
severely  restrict  the  range  of  possible  electron  density  profiles  which 
could  have  produced  those  measurements.  Our  previous  discussions  suggest 
that  the  density  profile  can  be  more  strongly  restricted  by  use  of  additi¬ 
onal  data  from  similar  lower  frequency  measurements  since  this  data  is 
influenced  disproportionately  by  the  peak  in  the  profile. 

In  order  to  carry  out  such  a  determination  of  the  ionospheric 
profile,  the  minimum  requirements  include  a  satellite  transmitting  very 
stable  plane  polarized  beacon  signals  at  several  frequencies  such  as  400 
MHz,  140  MHz  and  20  MHz  along  with  other  phase  related  signals  at  slight 
offsets  from  these  frequencies  for  differential  measurements.  (It  might  be 
necessary  to  have  a  wider  variety  of  frequencies  if  measurements  at  all 
times  are  to  be  useful  since  the  peak  density  shows  substantial  varia¬ 
tion.)  At  least  one  ground  receiving  station  is  needed  at  such  a  location 
that  the  signal  path  is  in  a  region  in  which  the  electron  density  is 
desired.  If  a  geostationary  satellite  is  not  used  the  interpretation  of 
the  results  will  be  more  difficult  because  horizontal  gradients  will 
influence  the  time  rate  of  change  of  u  or  However,  if  absolute  mea¬ 

surements  of  ii  and  can  be  obtained,  as  is  done  for  ATS-6,  the  hori¬ 
zontal  gradients  may  be  inferred  from  the  time  rate  of  change. 

Once  multiple  measurements  of  i2  and  i))^  are  available  at  sev¬ 
eral  frequencies,  and  since  these  values  depend  differently  upon  proper¬ 
ties  of  the  profile,  it  is  possible  to  match  the  experimental  values  by 
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variation  of  the  assumed  ionospheric  profile.  This  is  done  by  choosing  a 
profile,  by  ray  tracing,  and  at  least  for  the  lower  frequencies,  by  estab¬ 
lishing  the  n's  and  (ti^^'s  which  would  result  from  this  profile.  This 
procedure  is  then  repeated  with  different  profiles  until  a  match  with 
experimental  data  is  achieved.  It  is  very  important  that  the  model  pro¬ 
files  used  be  functions  of  a  relatively  few  parameters  which  influence 
different  features  of  the  profile  in  order  to  avoid  ambiguities  in  best 
fits.  These  parameters  may  be  strictly  geometrical  or,  as  with  Davies^^, 
a  mixture  of  geometrical  and  physical  parameters. 

The  works  of  Potts^^  and  Davies^’^  demonstrate  that  the  parameter 
upon  which  the  results  depend  most  strongly  is  the  peak  electron  density. 
Both  of  these  efforts  used  auxiliary  data  to  fix  this  peak  density.  This 
was  a  key  element  in  their  success.  Any  attempt  to  fix  an  electron  density 
profile  at  a  particular  time  and  place  will  probably  be  useless  unless 
there  is  an  experimental  method  to  establish,  with  confidence,  the  peak 
density  along  the  ray  paths.  Note  that  the  ray  paths  for  different 
frequencies  and  for  the  X  and  0  waves  are  different  so  that  horizontal 
gradients  must  be  taken  into  account  even  if  the  satellite  beacon  is  geo¬ 
stationary.  As  a  result,  some  means  would  be  required  to  estimate  these 
gradients.  This  could  be  done  by  using  rate  of  change  data  for  a  nonsta¬ 
tionary  satellite  to  find  gradients  along  the  satellite  track  or  by  using 
a  small  array  of  ground  receivers  spaced  on  the  same  scale  of  the  maximum 
deviation  of  the  ray  paths  from  straight  lines  which  may  be  tens  of  kilom¬ 
eters.  This  partial  redundancy  of  data  would  also  serve  to  discriminate 
against  atypical  cases  and  systematic  error.  Those  receivers  located 
along  a  line  perpendicular  to  the  moving  satellite  track  give  the  poten¬ 
tial  for  mapping  out  the  electron  density  profile  over  a  large  area  with  a 
single  pass. 

Methods  which  measure  phase  path  deviations  such  as  Faraday 
rotation  or  modulated  phase  are  subject  to  substantial  degradation  by 


87 


scintillations  imposed  by  irregularities  in  the  ionospheric  structure**^. 
Spread-F  and  sporadic-E  are  seen  to  produce  anomalous  phase  path  fluctua¬ 
tions  which  effectively  produce  an  increased  uncertainty  in  the  desired 
mean  value  of  such  quantities  as  a  and  (i)^.  As  a  result,  such  methods  as 
we  have  suggested  are  best  applied  only  during  times  when  the  ionosphere 
is  relatively  smooth.  This  is  especially  true  at  those  frequencies  near 
the  peak  plasma  frequency  since  scintillation  increases  strongly  with 
decreasing  frequency.  In  this  same  vein,  it  should  be  noted  that  the 
earth's  magnetic  field  is  a  dynamic  quantity'*^,  and  1  percent  variations 
in  Faraday  rotation  due  exclusively  to  field  fluctuations  are  typical. 
During  magnetic  storms  greater  variations  occur.  As  a  result,  the 
application  of  Faraday  rotation  data  to  ionosphere  profile  parameters  will 
always  yield  some  uncertainty.  As  indicated  by  Equation  6-15,  the 
measurement  of  the  path  integral  of  u  or  of  the  modulation  phase,  <!)„,,  is 
also  subject,  to  a  mild  extent,  to  second  order  variations  in  B  since  V  « 
B.  However,  it  is  possible  to  eliminate  the  first  order  Y  variations  by 
using  a  phase  deviation  (jin,  which  is  the  sum  of  the  left  and  right 
circularly  polarized  phase  deviations.  This  sum  includes,  to  lo»i«st  order 
in  Y,  a  term  proportion  to  Y^  which  should  be  negligible  since  Y  <  0.1  for 
even  the  lowest  frequencies. 
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SECTION  7 

CHEMILUMINESCENT  METHODS 


Atmospheric  ionic  and  neutral  species  can  be  produced  in  excited 
states  through  a  variety  of  processes  including  chemical  reactions,  inter¬ 
action  with  precipitating  electrons,  interaction  with  solar  photons, 
collisions  with  photoelectrons,  and  collisions  with  thermal  electrons  in 
the  high  energy  tail  of  the  distribution.  The  decay  of  these  excited 
states  produces  radiation  in  the  form  of  photons  of  energies  which  are 
specific  to  the  decaying  state.  Therefore,  the  rate  of  emission  of  such 
photons  is  directly  related  to  the  rate  at  which  the  excited  state  decays, 
that  is,  the  density  of  atoms  or  molecules  in  the  state  divided  by  the 
lifetime.  In  a  steady  state,  the  rate  at  which  the  states  decay  is  the 
same  as  the  rate  at  which  they  are  populated.  Typically  observations  of 
this  airglow  resulting  from  decay  of  the  excited  states  are  used  to  relate 
the  intensity  to  some  integrated  properties  of  the  radiating  region  relat¬ 
ed  to  the  rate  of  population  of  the  excited  state. 

While  airglow  methods  are  potentially  applicable  to  a  wide  vari¬ 
ety  of  atmospheric  regions  and  conditions,  it  is  our  primary  interest  to 
look  into  cases  which  can  yield  information  about  the  electron  density  of 
the  ionosphere. Thus,  the  altitude  regime  of  interest  is  from  100  km 
to  above  500  km. 

In  the  auroral  region,  the  major  source  of  excited  states  of  the 
constituents  is  the  precipitation  of  electrons  which  have  eV  or  greater 
energies.  These  electrons  serve  to  populate  excited  states  by  collisions. 
As  a  result,  airglow  techniques  in  auroral  regions  chiefly  provide  infor¬ 
mation  about  the  distribution  of  these  non-thermal  electrons. 
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Outside  the  auroral  oval,  in  daytime,  solar  radiation  and  the 
resulting  photoelectrons  constitute  the  major  mechanisms  for  population  of 
excited  states.  Therefore,  observations  of  daytime  airglow  (which  is  dif¬ 
ficult  due  to  the  high  background)  have  the  potential  for  giving  informa¬ 
tion  about  the  character  of  the  solar  radiation. 

Outside  the  auroral  oval,  at  nighttime,  chemical  reactions  are 
the  major  source  of  excited  states  which  radiate  in  the  visible  and  UV  so 
that  airglow  observations  can  be  related  to  species  densities  and  reaction 
rates  for  those  reactions  which  populate  the  excited  states.  Rate  con¬ 
stants  for  the  ionic  reactions  of  interest  have  been  studied  for  some  time 
both  in  laboratory  and  aeronomy  experiments^®  so  that  the  rates  are  fairly 
well  known.  Furthermore,  the  neutral  species  densities  are  much  less 
variable  than  the  ionic  constituents.  Models®®  exist  which,  given  suit¬ 
able  parameters,  can  provide  the  neutral  densities  perhaps  within  a  factor 
of  two.  As  a  consequence,  it  appears  possible  that  observations  of  the 
nightglow  radiations  excited  by  ionic  reactions  can  be  used  to  provide 
information  on  charged  species  densities.  Factors  involved  in  application 
of  this  chemiluminescence  to  evaluation  of  electron  densities  will  now  be 
considered. 

In  the  following  sections,  we  will  discuss  those  features  of 
atmospheric  chemiluminescence  which  can  be  applied  toward  evaluation  of 
ionospheric  electron  density  profiles  and  practical  techniques  which  will 
allow  this  to  be  done.  Sections  7.1,  7.2,  and  7.3  give  those  essential 
features  of  nighttime  chemistry  and  resultant  emissions  which  are  directly 
related  to  the  electron  density.  Atomic  oxygen  is  found  to  be  the  most 
useful  emitter.  Section  7.4  considers  the  restrictions  on  observations 
which  are  imposed  by  actual  emission  rates  while  Sections  7.5,  7.6,  and 
7.7  relate  to  the  capabilities  of  some  practical  satellite  sensing  tech¬ 
niques.  Finally,  the  prospects  for  application  of  these  methods  are  given 
in  Section  7.8,  and  a  summary  is  provided  in  Section  7.9. 


7.1 


NIGHTTIHE  F-REGION  IONIC  CHEMISTRY 


The  densities  of  the  prominent  neutral  species  in  the  F-region, 
0,  N2  and  O2,  are  shown  in  Figure  20.  The  corresponding  charged  species 
are  shown  in  Figure  8  to  be  O'*',  NO'*'  and  02^.  Chemical  reactions 
among  these  species  plus  electrons  have  the  potential  for  providing 
excited  electronic  states  of  any  of  the  six  major  constituents.  The 
reactions  which  will  be  considered^^"^^  are  given  in  Table  12.  This  set 
of  reactions  is  generally  thought  to  provide  the  significant  features 
which  fix  densities  of  the  ionic  species. 

All  of  the  reaction  products  can  potentially  be  created  in 
excited  states  and  thus  be  candidates  for  sources  of  emissions  which  can 
provide  information  about  the  constituent  densities.  However,  the  O2 
product  is  not  useful  for  establishing  ionic  densities,  because  it  will  be 
swamped  by  0+0  ♦  O2.  Atomic  nitrogen,  N,  has  an  excited  state  which  is 
produced  by  the  NO'*’  dissociative  recombination.  However,  the  lifetime  of 
this  state  is  26  hours  which  precludes  its  use.  Likewise,  the  charged 
species  are  not  promising.  O'*"  is  not  a  reaction  product  since  it  is 
primarily  formed  during  the  sunlit  hours  or  by  transport  in  conjunction 
with  charge  exchange  with  the  protonosphere.  Reactions  producing  0“, 
02’*’,  and  NO'*’  all  have  quite  limited  energy  available  for  production  of 
excited  states  so,  to  the  extent  to  which  such  states  exist,  their  popula¬ 
tion  by  thermal  electrons  in  the  high  energy  tail  of  the  distribution  will 
compete  with  that  by  ionic  chemistry.  There  appear  to  have  been  no  obser¬ 
vations  of  line  or  band  radiation  from  excited  states  of  these  three 
constituents  at  night.  We  are  thus  left  with  the  most  promising  source  of 
useful  nighttime  airglow  as  atomic  oxygen.  There  are  several  reactions  in 
Table  12  which  lead  to  0,  so  that  there  is  potential  for  populating 
excited  states  by  a  variety  of  chemical  routes.  This  variety  is  not 
desirable,  so  we  will  now  consider  in  detail  features  of  the  various 
reactions  to  find  excited  states  which  are  populated  by  a  subset  of  the 
reactions . 


neutral  species  densities  at  F-reglon  altitudes  for  exoatniospheric  temperatures 
I  and  700*K.  (From  Reference  54.) 


Table  12.  Ionic  reactions  important  in  nighttime  F-region  (adapted  from 
Reference  55). 


NO.  REACTION  RATE  (cm Vs) 


1. 

O2  +  e~  -»■  0  +  0  +  6.  96  eV 

+ 

O2  Dissociative  Recombination 

.  ,0-7  1  T  ,-.63 
2.1x10  - 

Soo 

2. 

NO"*^  +  e“  N  +  0  +  2. 76  eV 

NO*  Dissociative  Recombination 

.  .0-7  1  T  ,-0.4 

4x10  - 

Soo"^ 

m 

0^  +  O2  O2  +  0  +  1.53  eV 

+ 

O2  Formation 

2x10'^'  f-JLj-O'** 

■ 

300 

1 

0^  +  N2  +  N  +  1.09  oV 

NO*  Formation 

,  ^  ,  -12  ,  T  N-1 

1.2x10  - J 

Soo'^ 

5. 

0^  +  e  +  0  +  hv  (13. 6  eV  max) 

Radiative  Recombination 

3.5x10-^2  ^JLj-0.7 

6. 

0  +  e  +0  +  hv  (1.46  eV  max) 

Attachment 

1.3xl0'^^ 

7. 

0“  +  0'*'  -»■  0  +  0  +  12. 1  eV 

Ion-Ion  Recombination 

2.7x10-7  (_L_]-0.5 

300 

8. 

0  +O+O2  +  0  +3.7eV 

0  Associative  Detachment 

2xI0"^° 
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7.2 


POPULATING  ATOMIC  OXYGEN  EXCITED  STATES 


The  states  of  atomic  oxygen,  01,  are  given  in  Figure  21  and 
transitions  which  have  been  observed  as  atmospheric  emissions  are  indicat¬ 
ed  by  darker  lines.  Excited  states  populated  by  reactions  are  limited  by 
the  energy  available  as  well  as  by  selection  rules  which  depend  upon  the 
quantum  numbers  of  the  states. 

NO'*'  dissociative  recombination^^  has  sufficient  energy  to  popu¬ 
late  either  the  N(^D)  or  the  0(^0)  excited  states.  However,  direct 
excitation  of  0(^0)  is  forbidden  by  spin  conservation.  The  reason  this 
dissociative  recombination  is  included  is  that  it  appears  that  the  N(^D) 
atoms  quickly  collide  with  ground  state  0(^P)  to  form  0(^0)  so  that  the 
net  effect  is  the  production  of  0(^0).  This  state  has  a  lifetime  of  110 
seconds  and  decays  to  the  ground  state  by  emitting  the  red  6300A  line.  (A 
more  complete  consideration  of  NO'*'  dissociative  recombination  will  require 
noting  that  not  all  these  reactions  give  N(^D)  and  that  for  a  large  enough 
O2  density  N(^D)+02  NO-K)  serves  to  compete  with  0(  ^D)  production.) 

0?'*'  dissociative  recombination^^  has  an  excess  energy  of  nearly 
7  eV  so  it  can  leave  the  pair  of  oxygen  atoms  in  a  variety  of  states: 
(3p,3p),  (^P, ^D),  (3p, ^S),  (‘d,^D),  and  (‘D, ^S).  The  distribution  of  0 
among  the  three  states  is  somewhat  uncertain.  The  work  of  Zipf^^  suggests 
that  this  distribution  is  strongly  affected  by  the  vibrational  state  of 
the  original  02'*'.  It  appears  that  if  02'*'  is  vibrational ly  excited  to 
V  >  3  states  then  the  number  of  0(^S)  states  produced  per  recombination  is 
0.1.  This  is  in  agreement  with  atmospheric  experimental  results^^. 
The  number  of  0(  ^D)  states  produced  per  recombination,  kp,  is  also  sub¬ 
ject  to  uncertainty.  Coggen  and  Smith^**  give  l^5+kQ=  0.6  while  Chandra  et 
al^*^  used  k^+k^  =  0.4  and  Tinsley*’**  used  k^+k^  =  1.  The  sum  k^+kp  is  used 
because  of  the  rapid  cascade  of  the  0(*S)  state  to  0(*D)  with  a  lifetime 
of  0.74  sec.  It  is  this  transition  which  gives  the  5577A  green  line.  The 
0(*S)  can  also  decay  to  the  ground  state,  but  the  rate  is  low  so  that 
94  percent  of  the  *S  decays  lead  to  the  5577A  Mne. 
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0'*'  radiative  recombination  has  an  energy  excess  of  13.6  eV  so 
that,  energetically,  a  large  number  of  excited  states  may  be  populated. 
Tinsley  et  al.^^  calculated  the  rate  coefficients  for  the  most  common 
excited  states;  and  these  are  reproduced  in  Table  13.  The  authors  estimate 
that  the  rates  have  an  accuracy  of  10  to  20  percent  for  the  larger 
coefficients  although  a  factor  of  two  would  probably  be  more  realistic. 
The  rate  for  the  905-910A  transition,  0905.910,  is  for  the  continuum  to 
the  ground  state.  Transitions  to  singlet  states  are  spin  forbidden  so 
that  the  6300A  and  5677 A  lines  do  not  appear.  The  most  promising  lines 
from  radiative  recombination  will  be  1356A  and  7774A.The  F-region  is 
optically  thick  to  the  1304A  emission.  The  8447A  and  4368A  lines  are 
possible  but  they  are  relatively  weak. 

Table  13.  Partial  rate  coefficients  for  (T*'(**S°)  radiative  recombination 
yielding  emission  in  transition  listed  (units  are  10"^^  cm^ 
sec~^).  (From  Reference  53.) 


Transition  (A) 

700‘’K 

1000°K 

ZOOO^K 

905-910 

4.5 

3.8 

2.8 

1027 

1.1 

0.9 

0.5 

1304 

5.2 

4.0 

2.4 

1356 

10.6 

8.1 

4.8 

4368 

0.11 

0.08 

0.05 

7774 

7.6 

5.8 

3.4 

8447 

3.4 

2.6 

1.5 

Ion-ion  recombination  gives  an  excess  energy  of  12.1  eV  to  share 
between  the  two  product  oxygen  atoms.  Olson,  Peterson,  and  Moseley^^  have 
calculated  partial  rate  coefficients  for  this  reaction.  The  total  rate 
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coefficient  found  differs  from  the  experimental  value,  but  it  is  believed 
that  ratios  of  the  partial  rates  are  reasonable.  In  Table  14,  the 
fraction  of  the  total  rate  coefficient  is  given  based  on  Reference  61. 
These  results  are  generally  thought  to  be  within  a  factor  of  two  of  the 
correct  values.  However,  these  calculations  involve  approximations  in  the 
treatment  of  complex  atomic  processes,  and  the  detailed  effects  of  these 
approximations  are  unknown.  The  density  of  0"  available  for  ion-ion 
recombination  is  influenced  by  the  loss  of  0“  by  associative  detachment 
so  that  this  must  be  included  in  the  full  set  of  reactions  required  to 
establish  the  rate  of  production  of  excited  0  atoms. 


Table  14.  Fraction  of  0++  0“  ♦  0  *  0  reactions  which  lead  to  emissions 
listed,  [yy  =  2.7x10“^  /300/T  (cmVsec)]. 


Emission 

Fraction,  6 

1304A 

0.47 

1356A 

0.53 

7774A 

0.42 

6300A 

0.12 

8447A 

0.45 

EMISSION  RATES  OF  ATOMIC  OXYGEN  LINES  IN  THE  F-REGION 


In  order  to  get  a  feel  for  the  magnitude  and  variation  with 
altitude  of  the  01  emissions,  let  us  consider  typical  species  densities  as 
found  in  Figures  8  and  20.  A  set  of  densities  is  given  in  Table  15  for  an 
exospheric  temperature  of  800".  This  temperature  fixes  the  neutral  densi¬ 
ties.  The  electron  density  comes  from  Figure  8,  while  the  ion  densities 
are  adjusted  for  the  steady  state  in  such  a  way  that  the  total  ion  density 
equals  the  electron  density.  This  means  that 
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Table  15.  Sample  nighttime  species  densities  (ciir3).[a(b) 
indicates  axlO^]  (from  Reference  54). 


SPECIES  ALTITUDE 


_ 

150 

200 

250 

r - “ 

300 

350 

400 

450 

500 

B 

2.3 

(10) 

4.6 

(9) 

1.3 

(9) 

4.2 

(8) 

1.4 

(8) 

4.8 

(7) 

1.7 

(7) 

6.  1 

(6) 

B 

2.3 

(10) 

1.7 

(10) 

1.9 

(8) 

2.7 

(7) 

3.9 

(6) 

6.2 

(5) 

1.0 

(5) 

1.7 

(4) 

2.4 

(9) 

1.2 

(8) 

1.  1 

(7) 

1. ) 

(6) 

1.3 

(5) 

1.5 

(4) 

1.8 

(3) 
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[otj  “  1.2x10"**  [O2J  X  [0^]/[ej  (7-1) 


[NO'^J  =  1.7xlO"‘’[N2j  X  [0'*’J/lej 


(7-2) 


and 


1.3x10'^^  [0|  [ej 
2x10"^°  [0]  +  1.7x10"^  [o'*'] 


(7-3) 


Above  the  peak  of  the  F- layer,  the  steady  state  condition  as  well  as  the 
assumption  of  no  flow  of  constituents  will  probably  become  increasingly 
less  reliable.  Hence,  these  methods  cannot  be  used  to  provide  information 
above  500  km,  although,  in  principle,  it  is  possible  to  include  such 
effects  in  finding  a  dynamical  relation  between  the  various  densities  and 
rate  coefficients. 
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The  species  densities  can  now  be  used  to  provide  volume  emission 
rates  if  it  is  assumed  that  the  populations  of  the  various  levels  are  in  a 
steady  state.  In  this  case,  the  rate  of  emission  is  equal  to  the  rate  of 
population  of  the  levels.  The  emission  rates  in  photons/cm^-sec  are  thus; 


^1356  =  “1356  [0^]  [fi]  0*^3  Y7  [0  ]  ]  (7-4) 

E7774  =  07774  [0^]  [e]  +  0.42  Y7  [0^]  [0  ]  (7-5) 

E6300  =  0.12  Y7  [O^J  [0‘]  +  (k3+  kp)  Yi  [O2J  [ej 

+  1  X  Y2  [NO'^j  [e]  (7-6) 

E5577  ®  kj  X  Yi  [O2J  [e]  .  (7-7) 


(Here  the  ox  is  the  coefficient  for  x  A  emission  from  radiative  recombi¬ 
nation,  and  the  Yj  are  the  rate  coefficients  for  reaction  j  in  Table 
12.)  These  are  given  as  functions  of  altitude  in  Figure  22. 

If  the  atmospheric  and  ionospheric  densities  are  changed,  the 
emission  rates  will  change  accordingly.  At  a  steady  state 


E5577  *  [02] 

[0^] 

(7-8) 

E6300  “  IO2] 

[0^]  (  1  +  C  [N2j/[02]} 

• 

(7-9) 

When  O'*'  is  the  dominant  ion  (near  the  peak  of  the  F-layer),  these  two 
emissions  are  proportional  to  the  electron  density.  The  radiative  recom¬ 
bination  contributions  to  1356A  and  7774A  give 
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Figure  22.  Volume  emission  rates  of  several  01  lines  from  dissociative  recombination  (DR),  radiative 
recombination  (RR)  and  ion-ion  recombination  (I-I)  for  typical  F-region  parameters  given 
in  Table  13.  (Quenching  of  6300A  has  been  ignored  and  this  becomes  a  serious  omission 
below  200  km.) 


Err  •  le]  [0*] 


(7-10) 


This  is  approximately  proportional  to  the  square  of  the  electron  density 
so  that  emissions  are  greatly  enhanced  by  a  strong  F-region.  The  example 
given  in  Figure  8  of  a  nighttime  density  of  10^  cm“^  is  common  and 
represents  a  case  in  which  sensing  techniques  should  be  able  to  operate. 

7.4  EMISSION  INTENSITIES  REQUIRED  FOR  OBSERVATION 

In  an  optically  thin  volume  emitter,  the  intensity  of  radiation 
in  photons/cm^-sec  received  in  a  detector  of  small  solid  angle  Aft  is 

1=^  /  E(t)dt  (7-11) 

4ir  0 

where  E  is  the  volume  emission  rate  and  i  is  the  distance  from  the  detect¬ 
or  along  the  axis.  This  assumes  that  the  source  uniformly  fills  the  field 
of  view.  The  column  emission  rate 

J  =  /  E(ji)dJi  (7-12) 

0 

is  a  frequently  used  measure  of  the  source  strength  for  atmospheric  emis¬ 
sions,  and  it  is  given  in  units  of  Rayleighs  (R),  10^  photons/cm^-sec.  For 
the  examples  in  Figure  21,  the  6300A  emissions  give  «  100  R,  the  5577A 
line  gives  *•  lOR  while  those  emissions  provided  by  radiative  recombination 
yield  *«  0.1  R. 

The  ultimate  restriction  on  detection  of  atmospheric  emissions 
is  the  presence  of  sufficiently  large  numbers  of  photons  for  adequate 
statistics.  This  means  for  any  given  source,  for  a  10  percent  level  of 


101 


I 


statistical  uncertainty,  at  least  100  photons  are  needed  in  the  sampling 
interval.  At.  Thus,  this  confidence  level  requires 

100  <  J  •  —  •  At  •  A  (7-13) 

4it 


where  A  is  the  effective  area  of  the  detector.  This  means  that  there  is 
an  inverse  relationship  between  source  strength  and  spatial  or  temporal 
resolution.  As  an  example,  the  Atmosphere  Explorer-C  had  a  sensor  with  a 
field  of  view  of  a  3"  half  angle  cone  (Aft  =  8.6x10"^  sr) ,  a  detector  of 
0.5  an  diameter  (A  =  0.2  cm^),  and  a  sample  time  of  At  =  128  msec.^^ 
Therefore,  J  >  5.7  Rayleighs  would  be  required  for  a  ±  10  percent  level  of 
uncertainty. 

In  the  case  of  a  detector  which  is  mounted  on  a  satellite  and 
continuously  viewing  earthward,  both  field  of  view  and  sample  time  might 
be  increased  substantially  over  the  above  example  if  the  horizontal  varia¬ 
tion  of  the  source  is  not  too  great.  Consider  a  satellite  at  1000  km  cir¬ 
cular  orbit  (period  =  105  min)  with  a  detector  field  of  view  of  500x500  km 
at  F-region  heights  (Aft  ->  0.5  sr).  A  sample  interval  might  be  such  that 
the  subsatellite  point  moves  by  less  than  500  km,  say  100  km.  Then, 

At  •»  10  sec  will  maintain  the  spatial  resolution.  For  A  =  0.2cm^,  the 
minimum  column  anission  rate  which  can  be  detected  to  a  10  percent  level 
will  be  «  10“^  Rayleighs.  The  two  examples  given  serve  to  illustrate  the 
intensity  required  for  detection  and  resolution  of  atmospheric  emissions 
for  a  line  of  sight  which  is  vertical.  As  we  shall  see  in  a  later  sec¬ 
tion,  the  use  of  a  limb  mode  can  increase  the  intensities  substantially. 

7.5  USING  EMISSIONS  TO  FIND  PEAK  ELECTRON  DENSITY  AND  HEIGHT 

Chandra  et  al.^^  as  well  as  Tinsley  and  Bittencourt^®  have 
demonstrated  that  observations  of  both  1356A  and  6300A  emissions  during 
nigiictime  at  latitudes  below  the  auroral  oval  can  be  used  to  establish  the 
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peak  electron  density  of  the  F-layer  and  height  of  that  peak  if  a  shape 
for  the  profile  is  given.  Here,  we  will  outline  the  formalism  used  in 
Reference  60  although  that  in  Reference  59  is  essentially  the  same. 


For  a  vertical  viewing  direction,  let  us  take  the  electron  den¬ 
sity  profile  to  be  n^S(z)  where  n^  is  the  peak  electron  density  at  z  =  h^. 
Then  S(h^)  =  1.  The  column  emission  rate  of  the  1356A  line  is  then  given 
by 


2 

Jl356  =  % 


/  «use  S(z)d.  .  dz  (7-14, 


for  radiative  recombination  and  ion-ion  recombination  when  ng  <=  [O'*’ ] 
as  is  true  near  the  F-region  peak.  The  first  term  is  proportional  to 
n^2  while  the  second  term  is  important  only  when  it  is  proportional  to 
n^2.  That  is,  when  Y8[0]  ^  T?  ^be  second  integrand  is 

independent  of  n^  while  in  the  opposite  limit  the  second  integral  is 
small  compared  to  the  first.  As  a  result  J1356  is  very  nearly 
proportional  to  n^^  and  very  nearly  independent  of  any  other  species 
densities.  Therefore,  observations  of  0^35^  give  n^  when  the  profile 
S(z)  is  assumed. 

The  column  emission  rate  for  the  6300A  line  is  almost  exclusive¬ 
ly  from  dissociative  recombination  of  02''’  and  NO'*’.  The  state  of  01 
is  quenched  by  collisions  in  the  lower  F-region.  The  quenching  rate  is 
taken  as 

9(z)  =  \  [Nz]  +  Sg  [e]  (7-15) 

from  Nj  and  electron  collisions  where  =  7x10" ^ ^cm^/sec  and  Sg  = 
1.7x10"®  cm^/sec.  If  q  is  comparable  with  the  inverse  of  the  lifetime  of 
the  state  (t  =  110  sec),  then  deexcitation  of  these  states  is  impor¬ 

tant.  The  density  of  0'*’  can  be  related  to  [e]  in  a  steady  state  giving 
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(7-16) 


# 


where 


__jlJ _ 

1  +  B(z) 


B(z)  =  ^ 

TTi  [ej  Y2  Le] 


Therefore,  the  6300A  emission  rate  is 


(7-17) 


Jeioo  =  f  i(kstk„)  n„  /  Ik  (7-18) 


[1+B(z)]  [l+q(z)»T] 


where 


c 


T4 

(ks+kD)Tr3 


«>  0.03 


(7-19) 


and  fi  =  0.76  is  the  fraction  of  decays  which  lead  to  the  6300A  line. 
Jsaoo  is  thus  nearly  proportional  to  n,n  multiplied  by  a  function  of  the 
height  of  the  profile  peak  due  to  the  appearance  of  [O2]  in  the  integrand. 
Therefore,  the  ratio  J630o/(J use) is  chiefly  a  function  of  l\n  while 
i'iusej^'^^  is  chiefly  a  function  of  0^. 

Using  this  method  Chandra  et  al.^®  have  made  direct  comparisons 
with  other  data  on  electron  density  profiles,  and,  to  a  substantial 
extent,  they  are  in  agreement.  However,  there  are  several  sources  of 

uncertainty  in  this  approach  including  the  choice  of  atmospheric  models, 

departure  from  the  assumed  Chapman  profile,  and  the  accuracy  of  the 

chemical  rate  constants  involved.  In  spite  of  these  problems,  it  should 
be  possible  to  find  useful  information  on  subsatellite  electron  density 

profiles  at  night  below  the  auroral  oval,  especially  if  other 
instruments  are  used  for  calibration.  Such  analyses  are  not  presently 
carried  out  on  a  routine  basis. 
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7.6 


SCANNING  PHOTOMETERS  TO  FIX  ELECTRON  DENSITY  PROFILES 


Thomas  and  Donahue^ ^  have  used  the  vast  quantities  of  data 
gathered  by  an  orbiting  scanning  photometer  to  deduce  electron  density 
profiles  at  night,  below  the  auroral  oval,  by  a  deconvolution  of  many 
scans  of  the  same  region  from  different  positions.  Using  this  method, 
they  were  able  to  find  the  volume  emission  rates  in  the  plane  of  the  orbit 
of  the  060-6  satellite.  This  was  translated  into  electron  densities  in 
the  plane  of  the  orbit  by  methods  which  will  now  be  detailed. 

Emission  of  the  01  5577A  line  at  night  at  latitudes  below  the 
auroral  oval  results  almost  exclusively  from  production  of  the  'S  state  of 
atomic  oxygen  by  dissociative  recombination  of  02''‘  and  electrons.  The 
volume  emission  rate  in  a  steady  state  is 

E5577  =  ^21^5^3  [^2!  ]  (7-20) 

so  the  column  emission  rate  will  be 

J5577  ^2'<S  ^  ^3  [O2]  (7-21) 

where  the  integral  over  i  is  along  the  line  of  sight  of  the  detector,  and 
f2  =  0.94  is  the  fraction  of  decays  which  yield  the  55774  line.  Note 
that  no  correction  for  quenching  of  the  state  is  needed  due  to  its 
short  lifetime  of  0.74  seconds.  Figure  22  indicates  that  5577A  emissions 
are  small  below  150  kilometers.  However,  at  about  100  kilometers  another 
excitation  mechanism  becomes  available  to  oxygen  giving  rise  to  volume 
emission  rates  at  least  comparable  to  those  peak  values  in  Figure  22. 
This  alternate  excitation  mechanism  does  not  directly  relate  to  electron 
density,  so  it  is  of  no  interest  here.  This  means  that,  for  a  sight  path 
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which  goes  through  100  km  the  column  emission  rate  may  be  swamped  by  this 
additional  source.  This  will  make  it  very  difficult  to  deduce  information 
about  electron  densities  at  higher  altitudes.  Therefore,  observation  of 
5577A  emissions  will  be  most  useful  if  the  line  of  sight  does  not  extend 
down  to  the  E-region  of  the  ionosphere.  It  should  be  noted  that  5577 A 
observations,  just  as  for  those  at  6300A,  give  information  about  [O’*"  ] 
rather  than  [e],  so  one  must  either  assume  that  [O'*' ]  «  [e  ]  or  provide  a 
relation  between  the  two  as  was  done  by  Equation  7-16.  This  result  then 
gives 


‘^5577  =  i  ^3_p2]  [e.Ll^  (7.22) 

5577  2  S  J  1  +  B(£) 

The  OGO-6  satellite  used  a  scanning  photometer  to  view  the  5577A 
airglow  in  the  plane  of  its  polar  orbit  (apogee  1100  km,  perigee  400  km). 
The  scan  ranged  from  10°  to  25°  below  the  horizontal  plane  and  was  carried 
out  in  18.32  seconds.  The  field  of  view  was  4,2*  horizontally  but  only 
7,5  minutes  of  arc  in  the  vertical  direction  of  the  scan.  Each  18  second 
scan  consisted  of  127  positions  separated  by  7.5  minutes  of  arc.  There¬ 
fore,  for  a  satellite  altitude  of  500  km,  the  field  of  view  at  250  km  is 
about  85  km  x  2.5  km  midway  through  the  scan.  (This  gives  the  absolute 
upper  bound  on  the  possible  resolution  of  volume  emission  rates,  although 
we  will  shortly  see  that  the  practical  resolution  is  less.)  This  corres¬ 
ponds  to  a  solid  angle  of  An  =  1.6x  lO"**  sr  and  a  sample  time  of  At  = 
0.145  sec.  The  emitting  region  is  scanned  about  25  times  from  different 
positions  during  a  single  pass  across  the  subauroral  zone.  Thus,  a  total 
of  about  3000  data  points  are  taken  during  a  single  pass.  Thomas  and 
Donahue  divided  the  emitting  region  into  areas  2  km  in  altitude  by  0.5° 
magnetic  latitude  (-58  km)  horizontally  in  the  orbital  plane  and  at¬ 
tempted  to  find  the  volume  emission  rate  in  the  resulting  14,600  cells. 
(The  region  of  interest  was  about  200  to  400  km  in  altitude  and  ±  25*  in 
magnetic  latitude.)  Volume  emission  rates  in  the  cells  were  estimated  by 
beginning  at  the  highest  altitudes  giving  detectable  emissions  to  fix  the 
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rates  in  the  upper  cells.  Lower  scans  then  were  used  sequentially  to 
establish  the  rates  in  lower  altitude  cells.  The  resulting  emission  mat¬ 
rix  was  smoothed  and  repeatedly  rescaled  to  obtain  agreement  with  the 
measured  intensities.  The  resulting  matrix  is  not  unique  since  3,000  data 
points  are  used  to  define  14,600  volume  emission  rates.  Still,  the 
results  were  found  to  agree  with  measurements  to  within  10  percent,  and  it 
appears  that  the  main  features  of  the  emission  contours  (which  are  large 
compared  to  cell  size)  do  not  depend  upon  the  details  of  the  reduction 
procedure. 


This  technique  was  used  to  find  an  electron  density  profile 
which  could  be  compared  with  that  obtained  by  ground  based  radar  sound¬ 
ings.  The  results  indicate  that  v^ile  the  shapes  found  by  each  technique 
are  quite  similar,  the  amplitudes  differ  by  about  a  factor  of  two.  Thomas 
and  Donahue  produced  a  fit  by  altering  the  coefficient  of  Equation  7-22 
and  suggesting  that  the  rate  coefficient,  y3,  for  O'*’  +  O2  >  02'*’  +  0  should 
be  reduced  by  a  factor  of  about  two.  However,  it  is  not  clear  whether  Y3, 
ks,  or  [02]  uncertainties  are  most  important  in  producing  the  difference 
between  the  two  experimental  techniques  for  finding  the  electron  density 
profile.  Nevertheless,  it  is  very  encouraging  that,  given  some  calibra¬ 
tion  for  the  product  Y3ks  [O2J  by  separate  methods,  the  electron  density 
profile  in  the  orbital  plane  can  be  established  to  within  a  factor  of  two 
provided  that  the  intensity  of  the  emission,  or  equivalently  the  electron 
density,  is  sufficiently  great.  The  parameters  for  the  OGO-6  experiment 
made  observations  of  [e]  of  less  than  10^  cm“^  above  300  km  difficult 
while  values  as  low  as  10**  below  250  km  are  possible  although  difficult. 
As  a  result,  the  OGO-6  airglow  measurements  are  best  for  large  electron 
densities  and  a  low  F-region  peak.  A  modification  in  experimental  param¬ 
eters  would  be  required  if  the  full  range  of  nighttime  subauroral  iono¬ 
spheres  were  to  be  mapped.  Donahue®*’  has  suggested  that  a  scanning 
photometer  should  be  a  standard  part  of  experimental  satellites,  but  the 
recommendation  has  not  been  carried  out.  No  current  or  planned  and 
published  experiments  along  these  lines  now  exist. 


107 


7.7 


AN  ALTERNATE  SCANNING  METHOD 


A  sensor  with  a  limited  field  of  view  on  a  spinning  satellite 
will  scan  an  angular  profile  in  the  course  of  a  revolution.  Generally,  it 
is  not  possible  to  relate  the  angular  variation  of  radiance  directly  to 
volume  emission  rates.  However,  if  the  emission  is  from  a  spherically 
stratified  layer,  the  inversion  is  possible. For  a  column  emission  rate 
of 


J(r  )  =  2x10 
t 


-6 


I 


r  sin  z 
0  0 


E(r )dr 

sin2  Zq) 


1/2  ’ 


(7-23) 


where  r.  =  r  sinz  is  the  tangent  ray  height  for  the  sensor  at  radius  r 

too  3  J'  3  Q 

with  viewing  zenith  angle  z^,  the  volume  emission  rate  can  be  solved  for 
as 


E(r) 


19.1 1. 

It  dr 


00 


I 

r 


■^t 

r(r^  -  r^)  */  2 


(7-24) 


Therefore,  given  spherical  stratification  in  the  viewing  plane,  the  volume 
emission  rate  can  be  deduced.  This  requires  use  of  a  set  of  measurements 
of  the  column  emission  rates,  J,  with  sufficient  accuracy  and  resolution 
to  allow  the  inversion  in  Equation  7-24  to  be  carried  out.  Generally,  any 
experiment  of  this  nature  will  have  to  be  designed  to  take  into  account 
required  resolution,  spin  rates,  and  likely  levels  of  horizontal  gradents 
in  order  to  be  useful.  The  existence  of  such  horizontal  gradients  will 
introduce  errors  into  the  spin  scan  method,  but,  with  a  sufficiently  large 
number  of  scans,  it  may  be  possible  to  correct  results  to  provide  adequate 
emission  profiles.  Of  course,  once  emission  profiles  are  known,  it  is 
possible  to  convert  them  into  electron  density  profiles  using  atmospheric 
chemistry  to  find  the  rate  of  population  of  those  excited  states  whose 
emissions  are  being  observed. 


7.8 


PROSPECTS  FOR  FINDING  ELECTRON  DENSITY  PROFILES  BY 
CHEMILUMINESCENCE 


Past  experiments  indicate  that  measuranents  of  nighttime  airglow 
in  the  subauroral  region  can  be  used  to  find  electron  density  profiles. 
Presently,  no  satellites  exist  which  have  suitable  instrumentation  for 
this  purpose  but  there  is  nothing  particularly  unusual  or  difficult  about 
the  required  sensors.  We  will  now  describe  features  of  instruments  which 
would  be  suitable  for  this  task. 

7.8.1  Inversion  Using  Limb  Scanning  Photometers  in  Orbital  Plane 

In  order  to  hope  to  recover  an  adequately  resolved  set  of  volume 
emission  rates  from  a  series  of  column  emission  rate  measurements,  it  is 
necessary  to  have  a  sufficiently  small  field  of  view  to  resolve  features 
comparable  in  size  to  those  of  electron  density  variations.  This  typical¬ 
ly  means  vertical  resolution  of  »  10  km  and  horizontal  resolution  of  «  200 
km.  The  sample  time  must  be  long  enough  to  provide  good  statistics  (say, 
better  than  ±  10%)  but  short  enough  that  the  satellite  moves  a  distance 
less  than  the  required  horizontal  resolution.  In  the  case  of  5577A  emis¬ 
sions,  limb  column  emission  rates  of  as  low  as  J  •»  10  Rayleighs  are  avail¬ 
able  at  times  when  the  electron  densities  are  below  10^  cm"^.  If  the 

orbital  radius  of  the  satellite  is  too  great,  resolution  can  be  gained 

only  by  a  small  angular  field  of  view  so  it  is  best  if  an  orbit  near  the 
F-region  peak  is  used.  Using  a  circular  orbit  at  500  km  in  conjunction 
with  the  indicated  field  of  view  and  a  sensor  area  of  A  =  1  cm^,  the 
available  column  emission  rate  will  force  a  minimum  sample  time  of  1/4 
second  for  features  2000  km  away.  (This  distance  is  approx imately  that 
from  a  500  km  high  satellite  to  a  point  at  250  km  altitude  along  a  line 
tangent  to  a  spherical  shell  passing  through  this  point.)  Significant 
emissions  are  expected  from  altitudes  of  500  to  100  km  so  that  (500-100)/ 
10  =  40  sample  intervals  will  be  needed  to  scan  the  limb.  This  will 


require  ten  seconds.  In  ten  seconds,  this  satellite  will  move  •>  80  km 
which  is  less  than  the  required  horizontal  resolution.  Furthermore,  iono¬ 
spheric  changes  generally  occur  on  scales  of  many  minutes.  Therefore,  10 
seconds  is  an  acceptable  interval.  Such  a  scanning  photometer  will  obtain 
a  statistically  significant  number  of  counts  for  limb  column  emissions 
from  about  150  to  450  kin  for  both  the  5577A  and  6300A  lines.  Other  lines 
emitted  from  this  region  will  not  be  of  sufficient  intensity  to  provide 
useful  data. 

In  order  to  find  electron  density  profiles,  it  is  first  neces¬ 
sary  to  unfold  col  imn  emission  rates  to  give  volume  emission  rates.  This 
is  most  easily  done  in  the  orbit  plane  if  all  scans  occur  in  this  plane. 
As  shown  by  Thomas  and  Donahue^^,  multiple  scans  through  the  same  region 
from  different  satellite  positions  can  be  used  to  deconvolve  the  volume 
rates  from  col  unn  rates  and  thus  give  the  volume  emission  rates  in  the 
orbital  plane.  Our  example  of  10  seconds  for  each  scan  of  40  altitudes 
will  give  23,000  data  values  per  orbit.  With  this  data,  we  could  hope  to 
determine  up  to  23,000  volume  emission  rates  in  the  plane  of  the  orbit. 
This  corresponds  to  cells  10  km  in  altitude  by  about  80  km  in  horizontal 
dimension,  which  is  consistent  with  the  resolution  required  for  many  iono¬ 
spheric  structures. 

Such  an  inversion  is  essentially  that  used  by  Thomas  and  Donahue 
so  that  there  is  no  doubt  that  this  is  a  practical  method  of  analysis.  In 
order  to  convert  volume  emission  rates  to  electron  densities,  it  is  neces¬ 
sary  to  apply  dissociative  recombination  chemistry  to  relate  electron 
density  to  population  rate  of  those  states  which  decay  by  emitting  the 
5577A  and  6300A  radiation.  This  relation  involves  partial  rate  coeffi¬ 
cients  and  the  density  of  O2.  Both  of  these  are  subject  to  systematic 
uncertainties  which  could  be  resolved  by  constraining  the  solutions  for 
electron  densities  by  use  of  additional  data  from  other  sources  such  as 
radio  frequency  sounders.  Once  the  normalization  is  fixed,  the  conversion 


of  volume  emission  rates  to  electron  density  could  be  continued  with 
confidence. 


7.8.2  Limb  Scanning  Outside  the  Orbital  Plane 

Column  emission  rates  for  limb  scans  can  be  found  for  any  azi¬ 
muth,  not  just  in  the  orbital  plane.  However,  it  is  apparent  that  a  very 
large  number  of  such  scans  will  be  required  if  volume  emission  rates  over 
the  ionospheric  regions  visible  to  a  satellite  are  to  be  found.  It  would 
be  necessary  to  have  many  sight  paths  through  all  those  regions  where  the 
emission  rates  are  desired.  This  is  quite  impractical  for  those  distant 
points  at  the  edge  of  the  visible  region  perhaps  3000  km  away  from  the 
orbital  plane.  However,  if  scans  are  made  near  the  plane,  the  region  of 
deduced  volune  rates  may  be  expanded  due  to  the  continuity  of  atmospheric 
parameters.  Since  the  scale  for  horizontal  variation  is  “  200  km,  it  is 
expected  that  the  emission  rates  in  a  second  plane  parallel  to  and  about 
200  km  from  the  orbital  plane  can  be  found  if  a  second  series  of  scans  are 
carried  out  at  an  azimuth  near  the  orbital  direction.  A  suitable  azimuth 
would  be  e  »  200  km/2000  km  =  0.1  radians.  A  sequence  of  such  scans  at 
increasing  azimuth  angles  will  allow  the  addition  of  one  plane  removed 
another  200  km  for  each  azimuth  used.  Using  such  an  approach  would  allow 
the  volume  emission  rate  in  a  swath  parallel  to  the  satellite  track  to  be 
found.  Obviously,  the  data  handling  and  processing  would  be  cumbersome 
for  a  large  number  of  azimuths,  but  it  should  be  straightforward.  Since 
the  parameters  originally  chosen  for  a  scanning  photometer  were  based  on 
continuous  operation,  it  would  be  necessary  to  have  a  number  of  sensors 
equal  to  the  number  of  desired  azimuthal  positions  unless  the  characteris¬ 
tics  were  significantly  changed. 

A  second  approach  to  finding  volume  emission  rates  outside  the 
orbital  plane  can  be  used  if  horizontal  gradients  over  the  sampled  region 
are  small.  In  this  case,  as  shown  in  Section  7.7,  a  single  series  of  limb 
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scans  of  column  emission  rates  can  be  easily  Inverted  to  give  volume  emis¬ 
sion  rates  (which  are,  by  assumption,  functions  of  altitude  only).  The 
requirements  of  small  horizontal  gradients  is  difficult  to  satisfy  for  a 
satellite  in  500  km  orbit  scanning  ionospheric  heights.  In  this  case,  a 
typical  path  has  a  horizontal  extent  of  «  1000  km  along  which  contribu¬ 
tions  are  significant.  However,  as  indicated  earlier,  the  likely  horizon¬ 
tal  scale  of  structure  is  perhaps  200  km.  As  a  result,  routine  measure¬ 
ments  of  this  sort  are  not  likely  to  provide  an  accurate  measure  of  the 
distant  volume  emission  rates.  Still,  an  estimate  of  the  average  emission 
rates  over  the  sampled  region  can  be  obtained.  In  conjunction  with  other 
measurements,  this  may  be  applied  to  finding  an  estimate  of  the  electron 
density  profile  at  points  up  to  2000  km  from  the  satellite.  Of  course, 
the  conversion  (p  volume  emission  rates  to  electron  density  has  the  same 
uncertainty  associated  with  chemistry  rate  coefficients  and  neutral  densi¬ 
ties  as  does  the  case  for  measurements  in  the  orbital  plane. 

7.8.3  Electron  Densities  Beneath  the  Satellite 

As  discussed  in  Section  7.5,  a  single  observation  of  6300A  and 
1356A  (or  7774A)  emissions  can  be  used  to  provide  the  maximun  electron 
density  and  its  height  beneath  the  satellite  if  the  shape  of  the  profile 
is  known.  The  column  emission  rate  from  the  1356A  line  is  as  low  as  •> 
O.IR  for  observation  directly  toward  the  earth.  If  a  horizontal  field  of 
view  of  50  km  X  50  km  is  used,  then  a  solid  angle  of  Afi  •»  (50/250)^  = 
4x10"^  sr  will  provide  adequate  resolution.  An  effective  sensor  area  of  1 
cm^  will  then  prove  adequate  statistics  for  a  sample  time  of  less  than  one 
second.  The  6300A  line  is  sufficiently  intense  that  the  detector  just 
suggested  will  easily  be  adequate. 

The  required  assumption  of  profile  shape  is  quite  a  strong  limi¬ 
tation  and  a  source  of  uncertainty.  However,  if  this  method  were  used  in 
conjunction  with  a  photometer  scanning  in  the  orbital  plane,  the  two 
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methods  would  provide  complementary  Information.  It  is  expected  that  the 
scanning  photometer  would  provide  a  reasonably  accurate  shape  for  the  pro¬ 
file  while  the  uncertainties  of  the  method  would  be  largely  confined  to 
the  overall  normalization.  This  normalization  could  be  supplied  by  the 
earthward  looking  sensor  just  described  through  observation  of  6300A  and 
1356A  radiations. 

7.9  SUMMARY 

Observations  of  nighttime  subauroral  airglow  of  01  at  5577A, 
6300A  and  1356A  can  be  used  to  provide  routine  estimates  of  electron  den¬ 
sity  profiles  along  and  near  the  orbital  plane  of  a  satellite.  The 
required  instrumentation  is  quite  conventional.  A  limb  scanning  photome¬ 
ter,  operating  at  5577A  and  6300A,  with  a  narrow  field  of  view  gives  a 
series  of  column  emission  rates  which  can  be  inverted  to  provide  volume 
emission  rates.  These  volume  emission  rates  can  be  converted  to  electron 
densities  by  a  knowledge  of  the  dissociative  recombination  rate  coeffi¬ 
cients  which  serve  to  populate  the  states  of  atomic  oxygen  which  decay  to 
give  the  emissions.  This  conversion  to  electron  density  is  uncertain  to 
within  a  factor  of  '2  but  if  additional  measurements  are  available,  such 
as  from  ionosondes,  the  conversion  could  be  calibrated. 

These  limb  scans  are  most  useful  in  the  plane  of  the  orbit  since 
multiple  samples  of  all  regions  are  then  possible.  However,  additional 
scans  at  angles  somewhat  off  the  orbit  can  provide  the  electron  densities 
along  strips  parallel  to  the  orbit  once  the  orbital  plane  has  been  filled 
in.  The  continuity  of  electron  densities  allows  the  information  from 
these  scans  outside  the  plane  to  be  used  to  continue  the  density  estimates 
beyond  the  plane.  A  strip  several  hundred  kilometers  across  will  be 
accessible  to  practical  measurement  in  this  manner. 


The  scanning  method  can  be  augmented  by  a  fixed  sensor  which  is 
pointed  earthward  and  which  samples  emissions  at  6300A  and  1356A.  This 
data  can  be  used  to  provide  the  height  and  electron  density  at  the  F- 
region  peak  if  the  shape  of  the  profile  is  known.  It  is  expected  that  the 
scanning  method  will  provide  this  shape. 

Satellite  observation  of  airglow  from  atomic  oxygen  could  rou¬ 
tinely  provide  electron  densities  over  a  significant  portion  of  the 
globe.  While  this  is  applicable  only  at  night  in  the  subauroral  region, 
the  same  observations  can  be  used  to  provide  volume  emission  rates  without 
any  geographical  restriction.  Since  the  volume  emission  rates  are  worth¬ 
while  by  themselves,  such  an  experiment  should  be  of  interest  to  a  wide 
community  of  atmospheric  research. 


SECTION  8 
IN  SITU  METHODS 


Instead  of  positioning  a  satellite  far  above  most  of  the  iono¬ 
sphere,  it  is  possible  to  have  the  satellite  fly  through  the  peak  of  the 
F-layer  or  even  somewhat  below  the  peak.  Satellites  will  remain  in  a  200 
km  high  circular  orbit  for  a  significant  period  before  falling  to  earth 
and  have  a  many  year  lifetime  at  300  km.  This  latter  height  is  close  to 
the  peak  of  the  F-layer.  Flying  directly  through  that  portion  of  the 
ionosphere  of  greatest  interest  permits  direct  and  detailed  measurements 
of  many  ionospheric  parameters.  Not  only  can  electron  density  be  measured 
to  substantial  accuracy,  but  also  ion  density  and  composition,  plasma 
motion,  electron  and  ion  temperatures,  and  magnetic  field  variations. 
Table  4  lists  tne  in  situ  capabilities  of  the  ISIS  satellites.  These  fly 
well  above  the  peak  of  the  F-layer  but  still  within  sufficient  plasma  for 
in  situ  measurements  to  be  made.  The  Atmospheric  Explorer  satellites  fly 
much  closer  to  the  peak  of  the  F-layer  and  have  a  variety  of  in  situ  capa¬ 
bilities.  Figure  23  shows  a  sample  of  results  from  one  experiment  on 
board  Atmospheric  Explorer-E.  This  particular  record  shows  a  bite-out  in 
ion  density  occurring  as  the  satellite  passes  the  dip  equator.  Such  bite- 
outs  provided  the  first  conclusive  proof  that  bubbles  exist  in  the  equa¬ 
torial  ionosphere^ 

The  major  drawback  to  in  situ  measurements  is  that  they  are  only 
obtained  along  the  satellite  orbit.  To  obtain  anything  approaching  a 
three-dimensional  view  requires  many  satellites.  From  a  cost  standpoint 
alone,  this  is  not  feasible.  In  situ  measurements  are  valuable  when 
remote  techniques  are  inadequate,  when  high  accuracy  is  required,  or  when 
they  are  part  of  a  larger  ionospheric  spacecraft  such  as  ISIS. 
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Figure  23.  Atmospheric  Explorer— E  data  (from  Reference  65) 


SECTION  9 

BREMSSTRAHLUN6  METHODS 


In  recent  years,  it  has  become  possible  to  measure  the  low  enoroy 
x-ray  flux  caused  by  electrons  precipitatinq  into  the  atmosphere.  In 

some  cases,  it  has  been  possible  to  relate  this  "bremsstrahlung"  to 
measurements  of  the  precipitating  electron  flux,  as  well  as  to  optical 
observations  of  bright  auroral  arcs. Such  observations  are  of 
interest  here  because  the  precipitating  electron  flux  is,  under  some 

circumstances,  a  significant  source  of  ionization  in  the  ionosphere. 
Hence,  it  is  hoped  that  by  measuring  the  bremsstrahlung  in  a  particular 
area,  it  will  be  possible  to  compute  the  precipitating  electron  flux, 

which  in  turn  can  be  used  to  compute  the  ionospheric  electron  density. 

This  chapter  addresses  the  calculation  of  electron  flux  and  leaves  the 

calculation  of  electron  density  to  others. 

Figure  ?A  shows  the  observed  relationship  between  bright  auroral 
arcs,  bremsstrahlung,  and  precipitating  electron  flux.  Because  of  this 
strong  correlation,  it  is  also  anticipated  that  a  similar  correlation 
exists  between  the  diffuse  aurora,  bremsstrahlung,  and  precipitatinq  elec¬ 
tron  flux.  Diffuse  aurora  are  of  greater  interest  to  this  study  because 
they  cover  a  much  larger  area  of  the  high  latitude  regions  than  do  the 
bright  arcs  and,  therefore,  may  be  important  to  a  larger  region  of  the 
ionosphere.  Diffuse  aurora  involve  lower  electron  fluxes  and  lov/er  elec¬ 
tron  energies  than  do  the  bright  arcs. 
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Figure  24.  Simultaneous  DMSP  observations  of  bremsstrahlung  x-rays, 

auroral  arcs,  and  precipitating  electrons.  The  loi»er  charts 
show  a  cross  section  through  the  auroral  arc  at  the  top  (from 
Reference  3). 
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9.1  BACKGROUND  ASSUMPTIONS 

In  the  following  discussion,  we  consider  the  problem  of  determin¬ 
ing  both  the  primary  flux  and  spectral  energy  distribution  of  precipitat¬ 
ing  electrons,  on  the  basis  of  satellite  measurements  of  bremsstrahl ung 
produced  by  the  electrons  in  the  atmosphere. 

The  problem  may  be  outlined  as  follows.  Electrons  precipitating 
into  the  atmosphere  over  a  broad  area  can  be  characterized  by  a  function 
s(T),  el ectrons/(m^  sec  keV),  the  number  of  electrons  per  unit  area  per 
second  per  unit  energy  having  the  kinetic  energy,  T,  which  enter  the  top 
"surface"  of  the  atmosphere.  This  is  the  spectral  flux,  assumed  iso¬ 
tropic,  integrated  over  all  directions  in  the  downward  directed  hemis¬ 
phere.  This  electron  flux  produces  an  upward  directed  photon  flux 
spectrum  f(k)  photons/(m^  sec  keV)  at  the  top  of  the  atmosphere,  where  k 
is  the  photon  energy,  and  f(k)  is  the  integral  over  all  angles  in  the 
upward  directed  hemisphere.  A  satellite  mounted  detector,  essentially 
outside  the  atmosphere,  and  sensitive  to  all  photons  entering  from  the 
lower  hemisphere  would  measure  this  flux^  Broad  area  precipitation  and 
planar  geometry  of  the  atmosphere  are  assumed. 

The  relation  between  the  two  flux  spectra  can  be  expressed  as  an 

integral , 


f(k)  =  I  4.(k,T)s(T)dT  ,  (9-1) 

k 

where  all  the  complex  physics  of  electron  scattering,  bremsstrahl  ung  pro¬ 
duction,  scattering  and  absorption  of  bremsstrahl ung  photons,  etc.,  is 
contained  in  the  kernel  (|»(k,T).^‘’  Using  a  Monte  Carlo  method.  Reference 
74  has  calculated  il)(k,T)  for  the  conditions  assumed  above  in  the  energy 
ranges,  10  keV  <  k  <  1500  keV,  15  <  T  <  2000  keV.  The  behavior  of  this 
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function  is  shown  in  Figure  25.  Since  photons  are  only  of  equal  or  lesser 

energy  than  the  electron  which  produces  them,  «^(k,T)  vanishes  when  T<k. 

For  this  reason,  the  plot  uses  T-k  as  an  abscissa.  Further,  the 

extrapolation  of  all  curves  given  in  Figure  25  must  pass  through  the 

origin  on  this  plot.  The  special  feature  that  the  photon  flux  at  energy  k 

is  produced  entirely  by  the  part  of  the  electron  flux  at  energies  greater 

than  k  has  important  consequences  for  the  estimation  of  electron  fluxes 

from  photon  flux  measurements.  The  curves  of  Figure  25  do  not  extend  to 

low  enough  energies  in  either  k  or  T  to  be  of  maximum  utility  for  diffuse 

aurora.  Also  at  low  energies  (0.4  -  3.0  keV),  the  photon  spectrum 

contains,  in  addition  to  bremsstrahlung,  the  K  x-rays  of  nitrogen, 

a 

oxygen,  and  argon.  Calculations  of  auroral  bremsstrahlung  and  x-ray 
emission  at  photon  energies  of  0.1  to  10  keV  have  been  performed . 

These  papers  calculate  the  photon  spectrum  resulting  from  assumed  electron 
spectra,  not  the  more  general  kernel  4>(k,T)  of  Reference  74,  but  the 

methods  described  are  equally  useful  to  calculate  the  kernel  for  low 

values  of  k  and  T.  Therefore,  in  principle,  the  kernel  incorporating  the 
physics  of  the  process  is  known  or  computable.  Although  specification  of 
the  kernel  at  low  energies  would  require  new  calculations  to  be  made,  we 
assume,  for  discussion,  that  the  kernel  is  known.  The  essential 
difficulties  in  the  remote  sensing  problem  lie  elsewhere. 

9.2  THE  INVERSION  PROBLEM 

In  common  with  many  problems  of  remote  sensing,  the  estimation  of 
spectral  electron  flux  from  photon  measurements  can  be  cast  in  the  form  of 
a  Fredholm  integral  equation  of  the  first  kind.  This  is  the  form  of  Equa¬ 
tion  9-1  if  f,  and  ^  are  assumed  known  analytic  functions  and  s  is  the 
unknown  function  to  be  determined.  Except  in  restrictive  and  rather 
special  conditions,  which  fortunately  include  the  conventional  integral 
transforms  (Fourier,  Laplace,  Mellin,  etc.,),  a  unique  solution  may  not 
exist.  Moreover,  in  the  remote  sensing  problem,  the  function  f(k)  is  only 
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4.(k,T)  MeV 


known  at  a  finite  number  of  energies  k,  within  a  certain  measurement 

error.  In  these  circumstances  the  numerical  inversion  of  the  integral  of 
Equation  9-1  is  an  ill-posed  problem.  Many  functions  s(T)  containing  more 
and  more  rapid  oscillations  and/or  discontinuities  will  satisfy  the 

equation  for  a  given  finite  set  of  f's.  These  features,  and  several 

approaches  to  the  general  inversion  problem,  are  discussed  in  Reference 
77.  Successful  inversion  methods  all  impose  additional  conditions  upon 
s(T)  other  than  its  ability  to  produce  agreement  with  the  observations. 
These  additional  conditions  may  be  in  the  form  of  smoothness  and 

boundedness  requirements,  membership  of  the  function  s(T)  in  a  certain 
class  of  functions,  or  incorporation  of  theoretical  prior  knowledge  of  the 
form  of  s(T) . 

We  discuss  three  inversion  procedures.  A  finite  difference  or 
quadrature  approximation  and  a  variation  of  parameters  formulation,  which 
are  both  representative  of  conventional  approaches,  are  sketched. 
Finally,  a  procedure  based  upon  the  maximum  entropy  formalism  is 
examined.  This  latter  procedure  is  recommended  for  the  problem. 

9.2.1  Quadrature  Approximation 

This  method  consists  of  approximating  the  integral  equation 
(Equation  9-1),  by  its  discrete  analog,  a  matrix  multiplication 

m 

=  I  ’"''i )  s(T. )  w.  i=l,...,n  (9-2) 

1  j=l  ’  J  J  J 

or, 

m 

f,  =  J  "ijS,  (9-3) 

Where  the  measurements  are  tiie  n  values  of  fj  observed  at  energies  k.,-. 
The  quadrature  enters  in  the  determination  of  the  weights  wj  for  the 
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matrix  elements  M..=  d)(k.T.)  w.,  and  the  function  s(T.)  will  be  determined 
ij  1  J  J  J 

at  a  finite  number  of  points  T^.  ;  j=l,...,n;  m<n.  To  achieve  maximum 

definition  of  the  unknown  function  s(T),  take  m=n,  its  maximum  value. 

Further,  let  the  indices  i  and  j  run  over  the  same  set  of  energies,  i.e., 

k.  =  T.  when  i=i.  The  matrix  M..  is  then  triangular,  since  elements  with 
IT 

i>j,  or  are  zero.  Depending  upon  the  quadrature  or  weight 

function  one  wishes  to  adopt,  the  diagonal  elements  of  the  matrix  can  be 
made  either  zero  or  non-zero.  In  the  latter  case,  the  matrix  can  have  an 
inverse  and  n  values  of  Sj,  j=l,...,n  can  be  determined.  When  the 
diagonal  elements  are  zero  the  matrix  strictly  does  not  have  an  inverse, 
but  one  may  still  formally  calculate  n-1  values  of  s  j ,  j=l,...,n-l.  In 
either  case,  the  triangular  structure  of  the  matrix  permits  sequential 
calculation  of  the  s j ,  starting  at  the  largest  j  value.  Assume  diagonal 
elements  are  non-zero;  then. 


(M  )■!  f 
'  mm'  m 


(9-4) 


I 

l=k+l 


(9-5) 


Simple  and  direct  as  this  approach  may  be,  it  suffers  from 
several  disadvantages.  The  unknown  electron  flux  spectrum  can  only  be 
estimated  at  a  finite  number  of  energies  equal  to  the  energies  at  which 
bremsstrahlung  measurements  are  made,  which  may  imply  only  very  coarse 
resolution  of  the  spectrum.  An  attempt  to  increase  the  energy  resolution 
by  increasing  the  number  of  bremsstrahlung  measurements  increases  the 
dimensionality  of  the  matrix,  which  accentuates  the  effect  of  measurement 
errors  upon  the  inversion  calculation.  It  is  evident  that  the  sequential 
calculations  required  by  the  inversion  algorithm  rely  upon  an  increasing 
number  of  differences  between  measured  quantities,  as  the  calculation 
proceeds.  Finally,  there  is  a  large  inaccuracy,  and  some  outright  arbi¬ 
trariness,  in  the  replacement  of  an  integral  by  a  matrix  multiplication. 
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9.2.2  Variation  of  Parameters 


Given  a  set  of  photon  flux  measurements  f | ,  i=l,...,n,  the  varia¬ 
tion  of  parameters  method  consists  of  fitting  a  function  of  m  parameters, 
m<n,  to  achieve  the  best  least  squares  fit  to  the  measurements.  For 
example,  assume  the  electron  flux  spectrum  is  an  exponential  function, 

s(T)  =  Ne"“^  (9-6) 

of  two  parameters,  N  and  a.  This  function,  inserted  into  Equation  9-1, 
predicts  a  measurement  at  energy  of  f-j(N,a)  which  is  a  function  of 
the  two  parameters.  The  minimum  of  the  sum, 

F  =  "  [f.  -  f.{N.a)j2  (9-7) 

i=l  ^  ’ 

occurs  at  the  parameter  values  which  enable  the  chosen  function  to  best 
reproduce  (in  a  least  squares  sense)  the  observed  data. 

Unless  the  number  of  parameters  m  is  equal  to  the  number  of 
measurements  n,  it  is  impossible  in  general  to  fit  all  measurements.  How¬ 
ever,  there  are  measurement  errors  and  if  the  fi  are  within  the  measure¬ 
ment  errors  of  the  f ^ '  the  calculated  electron  flux  sperti^  jt  canno*^ 
improved  upon  within  the  family  of  functions  chosen.  This  method  yields 
the  solution  s(T)  most  compatible  with  the  observed  measurements  within 
the  family  of  functions  selected.  Its  success  depends  upon  judicious 
choice  of  the  functions  to  be  varied,  since  these  must  have  the  flexibil¬ 
ity  to  reflect  the  actual  structure  required  in  s(T).  If  one  accepts  this 
arbitrariness,  the  goodness  of  fit  to  the  data  can  be  quantified  by 
conventional  statistical  measures  applicable  to  least  squares  problems. 
However,  since  there  are  many  solutions  s(T)  which  will  reproduce  the 
observed  data,  goodness  of  fit  alone  is  no  justification  for  the  validity 


of  the  computed  electron  spectral  flux.  Unless  supplemented  by  other 
knowledge  (theoretical  reasons)  for  believing  the  electron  spectra  to  be 
exponential  for  example,  the  solution  obtained  depends  upon  the  function 
set  chosen  and  its  validity,  i.e.,  how  likely  it  is  to  be  correct,  cannot 
be  quantified. 

In  our  judgment,  the  limitations  of  both  conventional  inversion 
methods,  the  quadrature  approximation,  and  the  variation  of  parameters, 
are  obviated  by  the  approach  based  upon  the  maximum  entropy  formalism 
which  is  described  below. 

9.2.3  Maximum  Entropy  Formalism 

The  important  point  to  recognize  is  that  the  inversion  problem 
requires  the  estimation  of  a  continuous  function  (the  electron  spectral 
energy  flux)  on  the  basis  of  a  finite  number  of  measured  quantities 
(photon  fluxes  at  selected  energies),  which  by  the  form  of  Equation  9-1 
are  special  types  of  averages  over  the  desired  function.  Viewed  this  way, 
the  problem  is  identical  to  estimating  the  probability  distribution  of  a 
variable,  given  only  the  expectation  values  of  a  finite  number  of  func¬ 
tions  of  the  variable.  This  la^’er  problem  in  probability  theory  has  a 
solution  which  is  logical  and  rational,  but  not  unique.  Namely,  the 
distribution  with  the  maximim  informational  entropy,  as  defined  below, 
which  still  satisfies  the  constraints  imposed  by  the  given  expectation 
values,  is  overwhelmingly  the  most  probable  one,  and  the  optimum  inference 
to  be  made.  These  ideas  have  their  roots  in  information  theory.  They 
were  pioneered  by  Jaynes^®  and  have  been  developed  by  him  in  a  series  of 
subsequent  papers.  Reference  79  is  devoted  to  practical,  engineering 
applications  of  this  approach,  and  a  recent  review  is  provided  by  the 
conference  report  of  Reference  80. 

Maximum  entropy  methods  are  a  formalism  for  the  rational  use  of 
partial  information.  They  enable  problems  which  are  ill-posed  in  a  formal 


mathematical  sense,  in  that  the  given  information  does  not  determine  a 
unique  solution,  to  be  attacked  on  the  basis  of  inductive  logic.  The 
predictions  (or  solutions)  obtained  are  not  the  deductive  solution  of  a 
physics  problem  since  the  available  information  is  insufficient.  Rather 
the  predictions  are  the  result  of  a  form  of  statistical  inference,  which 
chooses  that  answer  which  has  maximum  statistical  weight  compatible  with 
the  physics  of  the  problem  (in  effect,  the  bremsstrahl ung  production 
kernel)  and  the  available  information  (the  measurements  of  bremsstrahl ung 
fl ux) . 

The  application  of  maximum  entropy  methods  to  the  bremsstrahl ung 
problem  and  the  inversion  of  Equation  9-1  is  a  new  use  of  the  formalism 
which  requires  some  special  analysis.  The  general  approach,  however,  is 
now  a  conventional  procedure  which  has  been  applied  to  problems  ranging 
from  estimation  of  power  spectra®^  to  optical  image  reconstruction.®^ 
Above  and  in  the  following  assertions  of  properties  of  these  methods  are 
made  without  proof.  The  requisite  demonstrations  can  be  found  in  the 
references  cited. 

9.3  THE  ESSENTIALS  OF  THE  MAXIMUM  ENTROPY  FORMALISM 

Following  Jaynes^®  we  explain  the  basic  idea  of  the  method  using 
a  discrete  probability  distribution.  Assume  that  a  quantity  x  can  take 
each  of  the  discrete  values  x^  (i=l,...,N)  with  a  probability  p-j .  In 
addition  to  the  knowledge  of  the  discrete  values  which  x  can  take,  there 
is  known  the  expectation  values  of  several  functions  fk(x),  k=l,...n, 
averaged  over  the  probability  distribution, 

\  =  I  JPi  k=l,...n  (9-8) 

Given  this  information,  the  probability  distribution  which  is 
least  biased,  in  the  sense  of  assuming  no  more  information  than  actually 


given,  is  the  distribution  which  maximizes  the  informational  entropy,  H, 
of  the  distribution  defined  as. 


H  =  -  5;  p.£n  p. 


(9-9) 


subject  to  the  constraints  of  Equation  9-8  and  normalization  of  the  proba¬ 
bility  distribution,  p-j  =  1.  The  constrained  maximization  is  effected 
by  associating  a  Lagrange  multipler  with  each  constraint  (Xq  for  normali¬ 
zation  and  for  each  of  the  expectation  values)  and  varying  the  p^. 
The  maximum  entropy  probability  distribution  is 

p.  =  exp[-Xo  -  [ 

The  Lagrange  multipliers  are  the  parameters  of  the  distribution. 
They  are  determined  by  the  corresponding  constraint  condition.  For 
example,  Xg  is  determined  by  the  normalization  constraint,  Xg  =  In  Z, 
where 


Z  =  I  exp[-5;  uJJx.)]  . 
i  k  ’ 

This  expression  is  recognizable  as  the  analog  of  the  partition  function  of 
statistical  mechanics. 

If  an  analytic  form  of  Z  is  known,  the  additional  n  equations  to  determine 
the  parameters  are 


fk  =  -(3/3Uk)ln  Z[ui,...,ij^) 


For  purposes  of  numerical  calculation  in  which  analytic  functions 
are  lacking,  the  procedure  is  to  directly  find  the  extremum  of  H. 


Equation  9-10  substituted  in  Equation  9-9  gives  H  as  a  function  of  the 
uk  s , 


=  In  +  1  u^f^  . 

The  values  of  uk  which  produce  an  extremum  of  H  are  the  correct 
parameters,  and  may  be  found  numerically.  Moreover,  the  function  H  can  be 
proved  to  possess  a  single,  unique  extremian  in  the  parameter  space. 

Fran  its  property  of  maximizing  the  entropy,  the  probability 
distribution  obtained  is  the  most  likely  distribution.  It  can  be  achieved 
in  more  ways  than  any  other  distribution  compatible  with  the  same  con¬ 
straints.  This  is  the  sense  in  which  maximum  entropy  methods  furnish  the 
most  logical  inference  from  the  information  available.  Whether  the 
inferred  distribution  is  or  is  not  correct,  is  a  matter  of  experimental 
measurement.  The  failure  of  direct  measurements  to  verify  a  predicted 
distribution  is  an  indication  only  that  the  initial  information  is  inade¬ 
quate— either  the  possible  values  |x,}  are  incorrectly  chosen  (essenti¬ 
ally  ignorance  of  the  possibilities),  or  the  functions  whose  average 
values  fk  are  given  are  not  sufficiently  sensitive  to  the  distribution. 
In  any  event,  the  maximum  entropy  prediction  is  the  best  estimate  possible 
on  the  basis  of  the  information  available. 

Despite  the  resemblance  between  Equation  9-8  and  the  bremsstrah- 
lung  production  Equation  9-1  there  are  certain  differences.  The  one  is  an 
integral  and  the  other  is  a  sum.  The  electron  flux  s(T)  is  not  a  proba¬ 
bility  but,  rather,  more  nearly  a  frequency  (number  of  electrons  of  a 
certain  type  in  the  sample).  The  next  section  derives  the  form  of  maximum 
entropy  inference  needed  for  bremsstrahl ung  applications.  The  resemblance 
of  this  exposition  to  standard  statistical  physics  is  not  accidental, 
since  the  same  methods  of  statistical  inference  are  common  to  both. 
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9.3.1  Derivation  of  Maximum  Entropy  Electron  Spectra 

To  employ  the  maximum  entropy  formalism  to  estimate  the  spectral 
flux  of  precipitating  electrons,  we  need  to  express  this  flux  in  terms  of 
a  normalized  continuous  probability  distribution.  We  define  the  single 
electron  energy  distribution,  p(T),  as  the  probability  per  keV  of  an  elec¬ 
tron  having  an  energy  in  the  range,  dT,  at  energy  T.  Then,  the  spectral 
flux  of  electrons  is, 

s(T)  =  G  p(T)  (9-11) 

where  G  is  the  total  flux  of  electrons  at  all  energies  (electrons/cm 
sec).  Bremsstrahlung  flux  measurements  at  n  different  energies,  k{keV), 
provide  n  equations,  derived  from  Equation  9-1. 

f|^  =  G  /  4,(k,T)p(T)dT  (n  equations)  .  (9-12) 

In  addition,  we  require  the  probability  to  be  normalized, 

1  =  /  p(T)dT  .  (9-13) 

We  maximize  the  entropy  of  the  single  particle  distribution,  since  we 
assume  electrons  are  independent.  The  appropriate  generalization  of  Equa¬ 
tion  9-9  for  continuous  distributions  is,®^ 

Hi  =  -  /  p(T)ln[p(T)/m(T)]dT  (9-14) 

where  m(T)  is  the  density  of  states  in  energy  for  precipitating  electrons. 

Because  of  the  presence  of  the  total  flux,  G,  an  additional 
unknown.  Equation  9-12  provides  only  n-1  constraints,  not  n  constraints. 
This  unknown  can  be  eliminated  by  taking  the  ratio  of  the  successive 
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equations,  indicated  by  Equation  9-12,  ordered  by  increasing  values  of 
observed  photon  energy,  k.  In  this  way.  Equation  9-12  is  re-cast  as  n-1 
constraint  equations, 

0  =  /  <t,(k+l,T)p(T)dr  -  ;  .^(k)p(T)dT  (9-15) 

where  k  and  k+1  indicate  successive  energies  at  which  bremsstrahlung 
fluxes  are  measured. 

We  now  maximize  the  entropy  given  by  Equation  9-14  subject  to  the 
normalization  constraint  of  Equation  9-13  with  Lagrange  multiplier  -a  and 
the  n-1  constraints  given  by  Equation  9-15  each  with  a  Lagrange  multiplier 
-3k,  k=l,...,n-l.  There  results  the  maximum  entropy  probability  distri¬ 
bution 


n-1 

p(T)  *  m(T)  exp(-a)  (9-16) 

The  factor,  exp(-a),  normalizes  the  probability  distribution.  Substitu¬ 
tion  of  Equation  9-16  into  Equation  9-14,  with  the  requirement  that  Equa¬ 
tion  9-15  be  satisfied  provides  an  expression  for  the  entropy. 


n-1 

Hi  =  In  /  dT  m(T)  [i'k'l’C^'^l  •  T)-Vi  ♦(k,T)]}  .  (9-17) 

The  values  of  6k  which  produce  an  extremum  of  Hi  are  the  correct  param¬ 
eters  to  use  in  the  distribution.  It  turns  out  in  this  case  that  Hi  is 
equal  to  a.  Once  the  distribution  p(T)  is  known,  any  one  of  the  original 
constraints,  given  by  Equation  9-12,  can  be  used  to  determine  the  remain¬ 
ing  unknown,  G. 
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There  remains  only  the  need  to  calculate  the  density  of  states, 
m(T).  The  number  of  single  particle  electron  states  in  which  an  electron 
has  a  direction  within  the  solid  angle,  dft,  and  magnitude  of  momentum 
within  dp  of  p  is, 

m  =  2ir  sine  de  p^dp  AXAyAz/h^ 

where  AxAyAz  is  the  volume  of  space  available  to  the  electron.  We  need 
the  number  of  states  available  to  an  electron  which  contributes  to  the 
downward  electron  flux.  All  electrons  within  a  distance  Az  =  jV^jAt  of 
the  "top  of  the  atmosphere"  and  with  angles  ii/2  <  e  <  n  contribute  to  the 
flux.  For  simplicity,  consider  non-relativistic  electrons,  so  that 
=  p  cos0/m  and  kinetic  energy  T  =  p^/2m.  The  integrated  number  of  single 
particle  electron  states  which  cross  area  AxAy  in  time  At  in  the  kinetic 
energy  interval  dT  at  energy  T  is, 

m(T)  =  (4xin/h^)T  dT  AxAyAt 

where  an  additional  factor  of  2  is  included  for  electron  spin  degeneracy. 
In  fact,  one  does  not  need  an  exact  expression  of  state  density.  Since 
both  p(T)  in  Equation  9-16  and  the  normalizing  factor  which  divides  it 
both  contain  m(T),  all  the  constants  divide  out.  We  need  only  know  that 
m(T)  is  proportional  to  T  in  order  to  calculate  the  maximum  entropy  esti¬ 
mate  of  spectral  electron  flux. 

9.3.2  Saaple  Calculation 

To  illustrate  the  foregoing  ideas,  we  present  a  simple  model  cal¬ 
culation  which  may  be  done  analytically,  but  which  preserves  the  general 
features  of  the  real  situation.  We  assume  a  simple  analytical  form  which 
roughly  approximates  the  true  kernel  of  bremsst rah lung  production  given  in 
Figure  25, 
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4)(k,T)  [keV-1]  =  3.703  x  lO'**  (T-k)  /  k^-’ 


T  >  k 


(9-17) 


=  0  T  <  k 

where  T  and  k  are  respectively  electron  and  photon  energy  in  keV.  We 
further  assume  an  electron  flux  spectrum,  s(T)  =  G  p(T),  given  by  the 
funct ion, 

p(T)[keV-i]  =  (T  -  bT2  +  cT3)  exp(-YT)  (9-18) 

with  the  normalization, 

N-^  =  2.273  X  lO-**  (keV)-^ 

and  constants, 

b  =  5.8  X  10*3  (keV)-^ 

c  =  10*®  (keV)*^ 

Y  =  10“^  (keV)*^ 

G  =  1.286  X  10^  elec/cm^-sec 

This  spectral  electron  flux  is  plotted  as  the  solid  line  in  Figure  26,  and 
designated  s(T). 

This  spectral  electron  flux  together  with  the  bremsstrahlung 
kernel  in  Equation  9-17  can  be  used  in  Equation  9-1  to  directly  calculate 
the  photon  flux  which  would  be  produced  at  various  energies.  The  results 
for  two  photon  energies  are: 

Photon  Energy,  k  (keV) 

50 
200 


1:2 


Flux  (photons/cm2-sec-keV) 
10.0 
0.49545 


We  now  reverse  the  logical  point  of  view  and  imagine  that  these 
two  photon  fluxes  are  measured  quantities  from  which  we  infer,  by  maximum 
entropy  methods,  the  electron  flux  h^ich  produced  them.  For  analytical 
simplicity,  photon  flux  measurements  are  considered  at  just  two  energies, 
although  the  more  measurements,  the  better  the  electron  flux  can  be 
determined.  With  two  measurements,  as  here,  only  two  parameters  of  the 
electron  distribution  can  be  established.  These  are  the  total  electron 
flux,  G',  and  a  single  Lagrange  multiplier,  g,  of  the  distribution  in 
Equation  9-16.  The  electron  distribution,  obtained  by  the  method  of  the 
previous  section,  consists  of  three  segments  because  the  exponential 
contains  zero,  one,  or  two  non-zero  terms.  (Remember  (|)  =  0,  if  T  <  k.) 
The  deduced  electron  flux  s'(T)  =  G'p'(T)  is  shown  as  the  dashed  line  in 
Figure  26,  which  may  be  directly  compared  with  the  original  electron  flux 
spectrum.  Analytically,  the  deduced  flux  spectrum  has  the  form, 

s{T')  =  G'p'(T)  =  G'A"^  T  exp{-g[f5o  4'{200,T)  -  f^oQ  (|,(50,T)]} 

(9-19) 


where  the  normalization  is, 

A-i  =  5.18  X  10-^  (keV)'2 


and, 

g  =  3.4235  X  10**  cm2-sec-(keV)  ^/photon 

G’  =  7.28  X  10**  (elec/cm^-sec) 

The  values  of  ^  are  given  by  Equation  9-17. 

In  this  example,  maximun  entropy  is  used  to  infer  an  original 
electron  spectrum  from  measurements  of  photon  fluxes  at  only  two  energies. 
This  is  asking  quite  a  bit.  But  Figure  26  shows  that  the  results  are 
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reasonable,  especially  at  higher  energies.  The  fact  that  the  prediction 
is  poorer  at  low  energies  is  not  surprising.  Because  the  breinsstrahlung 
production  kernel  is  zero  when  electron  energy  is  less  than  x-ray  energy, 
the  data  used  in  fact  supply  jio  information  about  the  electron  spectrum 
below  50  keV,  and  only  one  item  of  information  about  the  electron  spectrum 
between  50  and  200  keV.  Only  above  200  keV  is  the  electron  distribution 
sampled  by  both  measurements  of  photon  flux.  Although  this  example  used 
only  two  photon  energies  for  analytical  convenience,  no  difficulties  are, 
in  principle,  encountered  with  more  measurements,  and,  with  more 
measurements,  the  inference  becomes  more  accurate. 

To  implement  the  method  requires  more  accurate  knowledge  of  the 
production  kernel,  (ji,  especially  at  low  photon  energies  (10  keV  or  less) 
and  measurements  of  x-ray  flux  at  all  significant  energies,  including  the 
very  lowest. 

A  final  comment  can  be  made  about  the  effect  of  errors  in  the 
measurement  of  x-ray  fluxes.  In  terms  of  the  formalism  given,  the  effect 
of  errors  can  be  assessed  directly,  because  the  measured  fluxes,  f|^, 
enter  directly  into  the  probability  distribution  given  by  Equation  9-16. 
Consequently,  the  uncertainty  of  the  electron  spectrum  attributable  to 
measurement  error  can  be  directly  determined  by  computing  the  effect  of 
small  changes  in  the  f^. 
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